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Grant  SPC  00-4040 


Final  Report 

Studies  of  Intensive  X,  y  and  Optic  Emission  in  Runaway 

Breakdown  Process 

This  work  is  driven  by  intersection  of  two  lines  of  activities  in  studying  of  die 

electromagnetic  emission  generated  by  the  suprathermal  electrons.  We  have  addressed  the 

following  advanced  problems  in  the  grant  proposal: 

1 .  Theoretical  investigation  of  runaway  breakdown  process  in  stationary  electric  field. 

2.  Investigation  of  runaway  breakdown  phenomena  in  laboratory  installation  in  conditions 
of  cyclotron  resonance. 

3.  Investigation  of  an  influence  of  combined  effects  of  runaway  breakdown  and  cosmic 
rays  on  lightning  processes  in  thunderstorm  atmosphere. 

The  main  new  rezults  obtained  in  this  work  can  be  summarized  as  follows: 

1.  The  full  kinetic  theory  of  a  new  physical  phenomenon,  runaway  breakdown,  is 
developed.  The  theory  takes  into  account  both  ionization  and  bremsstrahlung  processes. 
The  kinetic  equation  is  fiilly  formulated  and  solved.  The  distribution  function  and 
avalanche  rate  are  determined.  The  electron-positron  pairs  generation  is  studied.  The 
possibility  to  observe  e^e“  aimihilation  line  during  intensive  high  altitude  discharges  in 
atmosphere  is  predicted. 

2.  Kinetic  theory  of  runaway  breakdown  in  inhomogeneous  thunderstorm  electric  field  is 
develped.  The  main  kinetic  equation,  boundary  and  initial  conditions  are  formulated.  Its 
solution  describing  the  electron  space  distribution  and  energy  spectrum  is  obtained. 
Spectmm  of  X-ray  emission  is  calculated.  It  has  quite  a  specific  form:  a  sharp  maximum 
at  energies  50  -  60  keV  and  a  rapid  fall  both  to  the  smaller  and  higher  energies.  These 
results  of  the  theory  are  in  a  good  agreement  with  observations.  The  intensity  of  the 
emission  is  growing  effectively  with  maximal  electric  field  Em  and  sharply  falling  down 
with  the  distance  from  the  point  z  where  E(z)=Em.  The  strong  growth  of  electrical 
conductivity  in  the  layer  AZ  near  this  point  is  predicted. 

3.  The  microwave  electron  cyclotron  resonance  discgarge  in  air  is  sstudied  in  a  magnetic 
mirror  trap.  The  obtained  experimental  data  demonstrate  a  rapid  growth  of  the  number 
of  fast  electrons  with  a  wide  energy  spectmm  (1-300  keV).  The  discharge  begins  with 
the  strong  X-ray  emission,  followed  by  the  syncrotron  emission  at  second  and  third 
cyclotron  harmonics.  The  optic  emission  from  a  low  energy  energy  part  of  electron 
spectmm  (<100  eV)  has  a  delay  up  to  80  ms  with  respect  to  the  X-ray  emission. 

4.  The  intensive  X-ray  emission  during  thunderstorms  was  studied  at  the  height  3340  m 
above  sea  level  by  the  specially  established  widely  spread  system  of  a  large  number 
(about  600)  SI5G  counters.  For  the  first  time  the  short-time  (1 — 5  min)  bursts  of  X-ray 
emission  were  observed.  The  bursts  are  highly  correlated  at  the  wide  space  region 
(about  0.5  km).  The  main  component  of  X-ray  emission  in  bursts  has  the  energy  50 — 
80  keV.  The  observed  bursts  could  be  attributed  to  the  bremsstrahlung  determinated  by 
the  mnaway  breakdown  (RB)  effect  in  the  thundercloud  electric  field.  Ground  based 
observations  of  RB  open  a  wide  opportunities  for  the  studies  of  fundamental  processes 
in  thunderstorms. 
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1 .  Theoretical  investigation  of  runaway  breakdown  process  in  stationary  electric 

field 


1. 1  Abstract 

The  full  kinetic  theory  of  runaway  breakdown  for  electrons  at  high  energy  range  is 
developed.  The  theory  takes  into  account  both  ionization  and  bremsstrahlung  processes.  The 
electron-positron  pairs  generation  is  studied.  The  possibility  to  observe  e^e“  annihilation 
line  during  intensive  discharges  in  atmosphere  is  predicted. 

1.2  Introduction 

Runaway  breakdown  is  a  new  phenomena  determined  by  the  effect  of  constant 
electric  field  at  high  energy  (s  >  100  Kev)  electrons  [1].  Its  physical  mechanism  is  deeply 
related  to  the  Coulomb  character  of  neutral  molecules  ionization  by  high  energy  electrons: 
collisional  crossection  a  in  nonrelativistic  region  drops  rapidly  with  energy  increase  a  oc 
7/s^ .  Thus  fast  electrons  in  the  presence  of  electric  field  E  >  Ec  {Ec  -  critical  field)  are  able 
to  become  runaways:  in  other  words  they  can  be  accelerated  by  the  electric  field. 
Accelerated  electrons  produce  a  number  of  low  energy  electrons  by  ionization  of  neutral 
molecules.  In  addition,  a  small  number  of  high  energy  electrons  may  be  generated.  These 
electrons  become  runaways  also  if  their  energy  exceeds  critical  value  Sc'-mc^ (Ec/E)  .  Being 
accelerated,  they  produce  another  generation  of  fast  runaway  electrons  and  the  process  is 
ongoing,  i.e.  as  a  result  an  avalanche  of  runaway  electrons  starts.  Exactly  this  process  is 
called  "runaway  breakdown".  It  is  determined  by  multiplication  and  acceleration  of  fast 
electrons.  The  exponential  growth  of  the  runaway  electrons  number  characterizes  the 
development  of  the  avalanche  and  absorption  of  electromagnetic  energy  by  discharge. 

The  kinetic  theory  of  runaway  breakdown  was  developed  previously  [2-5].  In  these 
papers  electrons  with  not  very  high  energies  {s<  5  MeV)  were  considered.  The  goal  of  the 
present  work  to  provide  the  full  kinetic  description  of  runaway  breakdown  effect  for  high 
energies  also.  We  will  demonstrate  that  at  high  energies  the  main  role  in  runaway 
breakdown  process  is  played  not  by  ionization  Coulomb  collisions,  but  by  much  more 
stronger  bremsstrahlung  effect.  The  structure  of  electron  distribution  function  at  high 
energies  would  be  determined.  The  role  of  energy  distribution  of  initial  fast  electrons 
triggering  the  runaway  avalanche  will  be  examined  in  details.  The  generation  of  electron  - 
positron  pairs  during  runaway  breakdown  will  be  studied  as  well.  It  would  be  demonstrated, 
that  the  last  process  is  quite  effective  what  allows  to  predict  the  possibility  to  observe  e'^e” 
annihilation  y-line  during  powerfull  atmospheric  discharges. 


1.3  Kinetic  theory  of  runaway  breakdown  at  high  energies 

At  high  energies  to  describe  the  runaway  breakdown  it  is  necessary  in  addition  to 
ionization  effect  take  into  account  bremsstrahlung  processes  for  electrons,  additional 
ionization  by  high  energy  photons  and  Compton  effect  for  photons.  To  estimate  the 
influence  of  all  these  processes  let  us  first  discuss  kinetic  equations  for  electrons  and 
photons  in  a  simple  "x"  approximation 
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Here  Ue  and  riph  -  number  densities  of  electrons  and  photons,  Nm  and  Z  -  number  density  of 
molecules  and  full  charge  of  their  nucleas,  a  =  e  /Pic  isa  fine  structure  constant,  ro=e  /me  - 
classical  radius  of  electron.  The  first  term  in  (1)  describes  the  absorbtion  of  photons  because 
of  Compton  effect,  4  is  Compton  absorption  length,  the  second  term  in  (1)  describes  the 
photon  generation  due  to  bremsstrahlung  (4^  is  the  length  of  generation  of  photon).  The  first 
term  in  (2)  describes  the  ionization  of  neutral  moleculas  by  fast  electrons  ( //  -  ionization 
length,  I(Z)=13.5  Z  eV),  and  the  second  one  by  photons. 

Let  us  first  compare  the  ionization  and  Compton  terms  (3).  Analyzing  the  equations  (1)  and 
(2)  it  is  possible  to  show,  that  in  conditions 


yln^(2y)  «  (n/aZ)  ln^[2'^^  mc^/I(Z)J 
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the  avalanche  is  determined  mostly  by  ionization  process.  In  this  work  condition  (4)  is 
supposed  to  be  fulfilled  always. 

The  full  kinetic  equation  describing  the  runaway  breakdown  process  takes 
inconditions 

(4)  the  following  form  [6]. 
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Here  f(p,Ht)  is  the  electron  distribution  function,  j?  -  electron  momentum,  ^  -  electric  field, 
p  =  (E  p}/(E  p)  -  cosine  of  pitch  angle.  The  integral  St  describes  collision  of  fast  electrons 
with  atomic  electrons  and  nucleus.  Fo  and  Fbr  is  fnction  forces  of  electron 
determining  energy  losses  due  to  ionization  and  bremsstrahlung: 


\dt !  me 

Here  Z,-  is  nucleus  charge. 

The  second  term  in  (5)  describes  pitch  angle  scattering  of  runaway  electrons.  Term  Si 
describes  appearance  of  a  new  bom  mnaway  electrons,  y  -  Lorentz  factor  of  electron,  8  and 
W- kinetic  energy  of  fast  electron  s  =  mc^(y-  1 ),  W=  mc^(y'-  1 ). 

At  Fig.  1.1  the  full  fnction  force  is  demonstrated.  One  can  see,  that  the  ionization 
plays  the  main  role  at  low  energies  s  <5  MeV.  At  £  >  5  MeV  the  dominant  role  belongs  to 
bremsstrahlung  force.  The  basic  relation  of  the  theory  which  determines  the  mnaway 
boundaries  eE  =  F  now  has  two  roots.  The  first  root  y  =  yi  determines  as  usual  critical 
electron  energy  Ec  [1].  The  mnaway  region  lays  between  yi  and  y2  .  The  second  root  is 
determined  mainly  by  bremsstrahlung  force.  The  mnaway  boundaries  energies  8c=mc^(yi-l) 
and  £c2=mc^(y2-l)  depend  on  relation  E/Ec 
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Fmin  is  the  minimal  drag  force  of  electrons  [1].  One  can  see  from  (7),  that  the  width  of 
mnaway  region 
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The  highest  energy  £c2  which  can  be  reached  by  mnaway  electron  is  growing  with 
E/Ec  approximately  linearly:  at  E/Ec  =5  it  reaches  the  value  430  MeV,  at  E/Ec=10  -  the 
value  850  MeV.  This  values  are  fiilly  determined  by  bremsstrahlung  force.  We  note,  that  in 
previous  papers  the  term  (6)  was  omitted.  It  means  that  kinetic  theory  was  considered  in  [1] 
for  low  energy  diapason  £«<£c2  only. 

Analyzing  collision  integral  we  see  that  pitch  angle  scattering  due  to  relativistic 
growth  of  electron  mass  with  energy  increase  falls  down  more  efficient  than  drag  force.  So, 
it  could  be  neglected  if  conditions 
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are  fulfilled.  At  such  conditions  pitch  angle  scattering  is  not  significant  and  due  to  this  the 
distribution  function  becomes  strongly  directed  along  electric  field  (see  [7]).  Thus  one  can 
consider  one  dimensional  -  integrated  over  angles  distribution  function. 
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-  bremsstrahlung  losses. 

The  ionization  term  takes  a  simple  form: 
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Thus  in  a  high  energy  range  due  to  relativistic  growth  of  the  electron  mass  its 
scattering  by  atomic  neucleus  is  weak  enough  and  the  distribution  function  of  energetic 
electrons  became  strongly  alongated  along  electric  field.  It  allows  to  use  kinetic  equation  in 
a  more  simple  form  (1 1),  (12). 


1.4  Time  development  of  distribution  function 

The  threshold  field  for  mnaway  breakdown  Ec  is  usually  much  less  than  for 
conventional  breakdown  Et.  As  an  example  for  air  and  for  hydogenium  Ec/Eth  ~  0.1.  On  the 
other  hand  to  excite  runaway  breakdown  one  needs  fast  electrons  with  energies  e  >  Sc.  That 
determines  an  important  difference  between  mnaway  and  conventional  breakdown.  How 
the  time  development  of  distribution  function  depends  on  the  existance  of  fast  electrons  or 
more  generally  on  initial  function  fo(yP  The  problem  was  studied  numerically.  The 
numerical  soutions  were  obtained  in  mnaway  region  (;'  <  yd)  by  predictor  -  corrector  and 
Crank  -  Nickolson  methods  with  boundaiy  condition  Both  methods  demonstrates 

identical  results.  To  determine  the  dependance  of  mnaway  process  on  the  energy  of  seed 
electrons  the  initial  function  fo(y)  was  taken  usually  close  to  delta  function  fo(y)°^6(y-y*). 

Results  of  numerical  calculations  of  equation  (1 1)  are  demonstrated  at  the  Fig.  1.2  - 
Fig.  1 .5.  We  see  from  Fig.  1 .2,  that  &iy*  <  yi  (or  e*  <  e)  independing  on  the  value  of  electric 
field  {E  >  Ec),  the  initial  distribution  function  dissapears  in  time  and  initial  number  density 
of  electrons  Ng  falls  down  exponentially.  So,  without  seed  fast  electrons  mnaway 
breakdown  does  not  exist. 

On  the  contrary,  if  y*  >  yj  the  mnaway  breakdown  is  effectively  developing.  Its 
main  features  seen  from  Fig.  1.3:  independently  on  the  value  y*  the  distribution  function  at 
low  energies  develops  very  fast  to  its  quasistationary  state 

f(yt)=Fo(y)e^'  (13) 

The  typical  time  of  this  process 

ti-X'  (14) 

At  the  high  energy  range  (y  >  y*)  the  establishing  of  quasistationary  solution  is  much  more 
slow  and  takes  time 


t2-dy2-y*)/(SE- 1) 


6e>1.2 


(15) 


Note  that  the  time  ta  is  diminishing  with  the  growth  of  y*  and  at  y*  >  ya  its  became  an  order 
of  ti. 

The  established  asympthotic  solution  takes  the  exponential  form  (13)  quite 
analogous  to  the  conventional  breakdown.  It  is  natural,  because  runaway  breakdown  is 
described  by  linear  equation.  Stationary  asympthotic  distribution  function  Fo(y)  is  shown  at 
Fig.  1.4.  for  different  values  of  parameter  Se-  One  can  see,  that  function  Fo  falls  down 
monotonically  with  y,  sharply  increasing  at  small  energies  e  =  mc^(y-l)  and  disappearing  at 
the  second  runaway  boundary  Fo(y2)  -  0.  The  behaviour  of  as5mipthotic  solution  Fq  could 
be  studied  analytically.  The  analytical  analysis  of  equation  (5)  shows  that  ionization  integral 
and  friction  term  are  most  essential  at  small  energies.  The  balance  of  these  terms  gives 
asymptotics  a.te-^0 
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In  intermediate  energy  range  yi  «  y  «  the  distribution  function  falls  down 
exponentially 
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To  determine  cutoff  the  distribution  function  near  right  stationary  point 

y2  we  can  neglect  the  ionization  term.  Expanding  friction  force  near  this  point  we  can  find 

the  solution 
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These  results  are  in  agreement  with  numerical  simulation. 


1.5  The  avalanche  rate 

One  of  the  most  important  parameter  of  breakdown  is  the  avalanche  rate.  The 
increasing  of  number  density  of  electrons  N(t)  is  shown  for  different  beam  energies  y*  at 
Fig.  1.5.  We  see  that  the  asympthotic  solution  demonstrates  exponential  growth  with 
parameter  X: 

N(t)=Noe^‘ 

We  investigated  the  dependence  of  avalanche  rate  X  on  dimensionless  electric  field  S 
.  The  results  are  shown  at  Fig.  1.6.  For  high  enough  values  of  S  this  behavior  is  well  fitted  by 
the  law: 

X  =  0.064  5^'^  (16) 

Note  that  as  seen  from  Fig.  1.5.  the  asympthotic  time  behaviour  of  electron  density 
does  not  depend  on  on  the  initial  distribution  function  (if  y*  >  y^  in  agreement  with  (13). 
Some  small  difference  in  li(t)  is  seen  in  the  begining  of  the  process  at  time  scale  of  the 
order  of  ti  (14).  We  see,  that  the  full  time  of  establishment  of  distribution  function  (15)  is 
not  significant  —  electron  density  growth  is  going  with  the  same  avalanche  rate  from  the 
very  begining.  It  means  that  the  avalanche  process  practically  does  not  depend  on  the 
established  distribution  function  Fo(y)  in  the  whole  breakdown  region.  The  dominant  role  is 
played  by  its  initial  part  in  the  vicinity  of  the  left  stationary  point  —  the  first  root  yi.  In 
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other  words,  only  the  region  near  critical  energy  8c  determines  the  avalanche  rate.  This  fact 
is  in  agreement  with  qualitative  theory  [1]. 

Basing  on  this  analysis  we  can  estimate  analytically  the  avalanche  rate.  Really,  near 
left  stationary  point  it  is  possible  to  consider  nonrelativistic  approximation.  In  this 
approximation  the  kinetic  equation  (5)  is  reduced  to  the  form: 


dt 


1  5/ 

-T  —  + 
du 


\fa{yv)wdw, 

u 


(17) 


There  is  an  important  property  of  the  equation  (17)  -  it  allower  to  eliminate  the  electric  field 
5e.  Introducing  variables 


T  =  2tS^,  v  =  u^J^ 

one  can  find 


%.+±± 
dr  dv 


dv  av^ 


\fo{yv)wdw 


V 


(18) 


Searching  exponential  solution  of  equation  (1 8)^0  =/ we  reduce  it  to  the  form: 


(1  _  ^  +  llo'  -  ^  =  0  (19) 

\a  )  Aay 

00 

Here  <I>  =  Jf(w)wdw,  and  variable  y  =  v^=  . 

u 

Note,  that  parameter  a  determined  by  (5)  could  be  cosidered  in  equation  (19)  as  a  large 
parameter.  For  a  Doo -we  have  from  (19) 

(l-y)O''-O'  =  0  (D  =  Coln(y-l)  +  C, 

This  solution  is  not  valid  for  y/Joo  and  yOl .  To  match  this  solution  let  us  return  to  the  finite 
big  values  of  parameter  a  and  use  an  expanding  variables  introducing  z  asy=  a^z. 

The  equation  (19)  takes  the  form: 

^ (1  -  -  {a^z  + 1)0,  -  =  0 

a  4a  z 

For  aOoo  its  solution  is 


-2az 


ho¬ 


using  a  procedure  of  matching  solution  [8]  we  find  "combined"  solution  which  is  good  with 
exponential  accuracy  everywhere  except  the  vicinity  of  pointy =7: 

(20) 


O  =■ 

y 


y-l 


In  the  vicinity  ofy=7,  introducing  x.’ y=  1+x  a  we  find  the  solution 

0  =  e 

It  means  that  solution  (20)  is  cut  at  u*^=l+  4/a.  As  we  shall  see  below  this  cutting  defines 
the  avalanche  rate  with  logarithmic  accuracy.  Integrating  (18)  from  u»  to  oo,  one  can  find: 


crln 


,16cr% 


4. 


1  2  L  4 

1H - h-  1-| - 

a  \  a 


+ 


1+-  +0^ 


(21)  8 


This  algebraic  equation  defines  parameter  X=2(y5^^/a.  For  a=ll  the  solution  of  (21)  gives 

X  «  0.055  5^^ 

This  result  is  in  a  good  agreement  with  numerical  calculations  (16).  More  than  that  it  is  in 
agreement  with  previous  numerical  solution  of  kinetic  equation  in  the  low  energy  limit  [2-5] 
and  with  approximate  qualitative  theory  [1].  The  reason  of  this  agreement  is  in  the  fact  that 
main  role  in  ranaway  breakdown  process  plays  the  low  energy  runaway  boundary  Sc- 


1.6  e*  e~  generation 

The  distribution  function  of  runaway  electrons  falls  down  rather  slow  (see(20)).  It 
means  that  fast  electrons  generates  intensive  flux  of  X  and  y  photons  due  to  bremsstrahlung 
process.  The  high  energy  photons  can  generate  electron  -  positron  pairs. 

The  process  of  X  and  y  generation  by  fast  electrons  in  neutral  gas  is  well  known. 
The  total  number  of  photons  could  be  estimated  using  emissivity  qy  [9] 


dN^  q^ 

-^  =  ^cNN, 
dt  s'”" 


■  AZ  co-f^e 


Here  Ng  -  the  number  density  of  fast  electrons  which  generate  bremsstrahlung  emission  and 
Nm  -  density  of  neutrals. 

The  crossection  of  e*  generation  is  quite  analogous  to  Bethe  and  Heitler 
bremsstrahlung  crossection  [9] 


da  =  ^ . smO^dO^, s{nG_d9_d<l)y.D 


In  coi^q* 


(22) 


K^K 


-{p ^  sin  9^  +  pi  sin  9_)- 


^P^P- 


K^K 


{2sl  +  2slq  )s\n9^.  sin^_  cos^ 


Here 

=£±-p±cos^±,  ={p^+P--kf, 

Integrating  this  crossection  we  can  estimate  the  e*  production 


f ++£•_=  6) 


dN 


dt 


^  =  a^cN^N^ 


a-.. 
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Taking  into  account  that  grows  exponentially  with  t,  we  find : 

lUz 

9  N, 


(23) 


This  relation  determines  the  number  of  generated  e^e”  pairs. 

To  estimate  the  emission  intensity  of  the  e^e“  annihilation  line  intensity  generated 
during  thunderstorm  it  is  natural  to  compare  with  the  results  of  observations  y  -  burst  from 
high  altitude  lightning  at  "Compton"  observatory.  According  to  Fishman  et  al  [10]  and 
Nemiroff  et  al.  [11]  the  flux  of  gamma  quants  in  energy  diapason  300  1000  KeV  observed 

by  BATSE  at  the  heights  500  km  was  of  the  order  of  100  photons/cm^s.  Using  relation  (23) 
and  taking  into  account  that  N/Ng  ~  a~  l(f^  one  can  estimate  the  flux  of  photons  from  e’^e” 
annihilation  at  the  distance  100  km  from  the  source  as  20  photons/cm^s.  Taking  into 
account  that  the  width  of  unihilation  line  is  narrow  one  can  suppose  that  the  y-quants 
generated  during  thunderstorm  could  be  detected  in  experiments. 
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Figure  1.1  Dependence  of  the  full  friction  force  F  on  y  (solid  curve).  The  dashed  curev 
presents  the  friction  force  without  accounting  bremsstrahlung. 
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2.  Kinetic  theory  of  runaway  breakdown  in  inhomogeneous  thundercioud  eiectric 

field 


2. 1 1ntroduction 


Multiple  observations  at  balloons  and  airplanes  of  electric  field  and  X-ray  emission  inside 
thunderclouds  revealed  a  definite  agreement  with  the  predictions  of  runaway  breakdown 
(RB)  theory  [1]: 

1.  Maximal  values  of  observed  in  thunderclouds  electric  field  \E\m  are  close  to  RB  critical 
field  [1,2], 


N„{z)  kV 


200^=4-4  — 


(1) 


nM  m 

Here  Nm(0)  the  density  of  molecules  at  the  sea  level  (Nm(0)=2. 7x  10^^  crn\  and  Nm(z)  -  the 
same  density  at  the  given  height  z. 

2.  The  observed  values  of  \Em\  fall  down  exponentially  with  the  height  z  in  accordance  with 
formulae  (1)  [2,3]. 

3.  The  definite  correlation  of  lightning  activity  with  the  regions,  where  electric  field  |.E|  is 
close  to  Ec  is  seen  [3]. 

4.  At  the  same  region  X-ray  bursts  -  strong  increase  of  X-ray  emission  during  one  minute 
time  intervals  are  observed  [4-6]. 

Note  that  the  structure  of  electric  fields  in  thunderstorm  atmosphere  is  strongly 
inhomogeneous.  The  thunderclouds  have  a  layered  structures  [3]  what  means,  that  the  main 
field  inhomogeniety  is  always  along  the  height  z 


E=E(z) 


It  is  of  importance,  that  not  only  the  value  of  electric  field  \E(z)\~Ec,  but  the  scale  of  field 
ingomogeniety  L  affect  strongly  the  RB  process  -  L  should  be  much  larger  than  the 
runaway  avalanche  length  U 

L»l„  (2) 

la  is  of  the  order  70  -  100  m  inside  thunderclouds  and  condition  (2)  is  usually  fulfilled  there. 
The  goal  of  the  present  paper  is  to  develop  a  kinetic  theory  of  RB  in  inhomogeneous 
thunderstorm  electric  field.  The  main  kinetic  equation,  boundary  and  initial  conditions  are 
formulated  in  section  2.  Its  solution,  presented  in  section  3  describes  the  electron  space 
distribution  and  energy  spectrum.  Spectrum  of  X-ray  emission  is  calculated  in  section  4.  It 
has  quite  a  specific  form:  a  sharp  maximum  at  energies  50  -  60  keV  and  a  rapid  fall  both  to 
the  smaller  and  higher  energies.  The  intensity  of  the  emission  is  growing  effectively  with 
parameter  EJEc  and  falling  down  sharply  with  the  distance  near  the  point  z  where  E(z)=Ec. 
Note  that  the  role  of  thunderstorm  cloud  inhomogeniety  was  discussed  in  a  simplified  form 
previously  in  [7-9] 

2.2  Kinetic  equation  in  inhomogeneous  electric  field 

Fast  electrons  generated  during  RB  process  in  nonuniform  electric  field  at  stationary 
conditions  (d/dt  =  0)  are  described  by  kinetic  equation  [10]: 


14' 


5//J  v"  dv 


dz  m 


v^dMV  m 

Here  f=f(v,nz)  -  is  the  electron  distribution  function,  //  -  is  a  cosine  of  an  angle  between 
direction  of  electric  field  E  and  velocity  v,  Fd  -  fiiction  force  due  to  ionization, 

V  ym) 

A  -  Bethe's  logarithm,  Vn  is  angle  frequency  scattering  parameter 


7^ 

^=Y— 

t!'4Z 


where  effective  charge  of  molecules,  Z,-  -  the  charge  of  every  moleculae  component: 

z  =  E,z.- 

The  term  Si  is  ionization  integral  [11].  Concrete  form  of  Si  will  be  presented  later. 

Electric  field  in  a  real  atmospheric  conditions  increase  inside  thunderstorm  cloud  up  to 
maximum  value  Em  and  then  decrease.  In  our  calculations  we  choose  it  in  a  parabolic  form: 


-L<z<L 


Here  Z,  is  a  scale  of  thunderstorm  electric  field. 

According  to  the  theory  of  runaway  breakdown  process  the  distribution  function  of  fast 
electrons  could  be  considered  approximately  as  strongly  directed  along  electric  field  [10]. 
Integrating  (3)  over  pitch  angles  d/j.  we  find  a  kinetic  equation  for  distribution  function  f(Y,z) 
of  fast  electrons 

(r'  (4) 

where  0d  is  normalized  relativistic  fiiction  force  [10],  E(z)=5(l-a^z^)  -  electric  field 
normalized  on  critical  electric  field  5  =  EJEc.  Height  z  -  normalized  on  ionization  length  4 
and  parameter  ai=la/L,  qcr(Y)  is  the  source  of  cosmic  ray  secondaries  generated  at  the  given 
height  z  in  atmosphere  by  the  flux  of  cosmic  rays  [12]. 

Ionization  integral  according  to  [10,1 1]  is  determined  by  formulae: 

(5) 

To  solve  the  equation  (5),  (6)  it  is  necessary  to  formulate  initial  and  boundary  conditions  for 
distribution  function.  As  is  well  known  runaway  breakdown  take  place  in  the  presence  of 
seed  energetic  electrons  only.  In  the  thunderclouds  the  role  of  seed  fast  electrons  is  played 
by  the  cosmic  ray  secondaries.  So,  the  initial  function  fo(Y)  is  the  energy  distribution  of 
cosmic  rays  secondaries.  It  is  determined  by  solution  of  equation  (\ref{fgam})  in  the 
absence  of  the  z  -  transport  term  {j^/dz)  and  electric  field  E(z).  The  same  function  fo(Y) 
determines  the  boundary  condition  at  the  initial  point  z=ai''.  The  natural  boundary  condition 
in  energy  space  is  f(Y)->0,  if  y^oo. 

Below  we  will  describe  the  results  of  solution  of  equation  (5),  (6)  for  the  given  initial  and 
boundary  conditions. 
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2.3  Distribution  function  of  electrons 

The  solution  of  equation  (5), (6)  was  studied  analytically  and  numerically.  The  results  are 
presented  at  Fig.2.1  -Fig.2.4.  At  Fig.  2.1  the  dependence  of  distribution  function  f(y)  on  the 
height  z  is  demonstrated.  One  can  see  that  the  growth  of  electric  field  E(z)  is  accompanied 
by  the  significant  growth  of  distribution  functionyf)^.  It  reaches  maximal  value  atz=2.5  and 
after  that /begin  to  diminish  with  z.  The  dependence  F  on  y&t  j<2  is  quite  well  described  by 
analytical  formulas  [10].  The  dependence  of  distribution  function  on  maximum  electric  field 
Em  is  shown  at  Fig.  2.2.  The  avalanche  type  growth  of  f  and  acceleration  of  electrons  at 
Em/Ec>l  is  seen.  The  same  avalanche  process  is  demonstrated  at  Fig.2.3  for  the  number 
density  N  of  runaway  electrons.  Note  that  at  £'=0  the  number  density  No  describes  the 
background  cosmic  ray  secondaries.  To  eliminate  this  background  which  is  not  connected 
with  RB  process  we  normalized  N  on  Nq.  At  Em>Ec  an  avalanche  growth  of  N  with  z  is 
clearly  seen.  Its  sharp  cut  off  coincides  with  the  points  z  where  E(z)  becomes  less  than  Ec. 
The  growth  of  the  maximal  avalanche  value  Nmax  on  Em/Ec  shown  at  Fig.  2.4  is  of  a  sharp 
exponential  form.  In  the  same  time  the  point  Zmax,  where  N(z)  reaches  its  maximum  value, 
moves  in  z  direction  with  Em/Ec  rather  smoothly,  reaching  for  example  z=2.6  4  at 
Em/Ec  =1.2  (Fig.  2.4). 


2.4  X-ray  emission:  energy  spectrum  and  space  distribution 

High  energy  electrons  passing  through  the  neutral  gas  generate  X-rays  due  to 
bremsstrahlung  process.  The  bremsstrahlung  intensity  and  spectra  is  determined  by  electron 
distribution  function/^f;  z): 


Jv  =  \<^v{s)wf{s,z)ds  =  ct.NJ^ 


(7) 


where  v  is  the  velocity  of  electrons,  a/e)  the  cross-section  of  the  production  of  photons  in 

2  2  2  2 

the  frequency  range  dv  by  incident  electron  of  energy  e  and  o2,=16/3-Z  a  ao  ,  ao=e  /me  . 
Using  expression  for  cross-section  [13],  one  can  find  flux  of  photons  jv: 


-1  +  -1-2v7^  ^ 

-\-2vy\ 


f{y,z)dy 


(8) 


The  expression  (7)  describes  the  source  of  bremsstrahlung  photons  at  given  point  z.  The 
dependence  of  the  flux yVon  the  energy  of  X  -  ray  quants  Ev  for  distribution  function/);  z) 
at  Em/Ec=l.b5  is  shown  on  Fig.  2.5  for  different  space  positions  z.  One  can  see  from  the 
figure  that  flux yV  as  a  function  z  is  initially  growing  with  z,  but  at  z=l  .375  reaches  maximal 
value  and  then  falls  down.  All  the  changes  of  yV  are  fully  determined  by  distribution 
function  of  runaway  electrons  and  its  dependence  on  the  electric  field  E(z). 

The  bremsstrahlung  photons  propagating  through  the  neutral  gas  loose  their  energy  due  to 
Compton  effect  and  disappear  due  to  photoionization.  The  kinetic  equation  for  photon 
distribution  function  I(v,z),  describing  these  processes  could  be  presented  in  a  form: 


clj  d^I  c  d 
Ij.  dv 


{v^l)-ca^J  +  J^=0 


(9) 


III  dz^ 

Here  z  is  as  previously  normalized  on  the  ionization  length  4  determined  by  RB  process, 
/r=7/2VeCTr  is  Thompson  scattering  length,  and  Oph- photoionization  cross-section. 

The  solution  of  diffusion  equation  (9)  with  the  source  (7),  (8)  is 

Vov(z 


l{v,z)-- 


h 

?yf2^ 


Vo-v 


■expi 


4D{vq  -  v)  9cr 


(10) 


Here  D=  Ij/Sla  and  we  take  into  account  that  for  a  considered  energy  range  £v>  20keV  the 
photoionization  crossection  could  be  chosen  in  a  form  [15]: 
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cr  .  =2*— OfloZ’ 

ph  3  0 


\mc  J 


^  7/ 


=  (JqV 


Here  as  usually  a=  e^/c  =  1/137  and  ao=e^/mc^. 

The  solution  (10)  is  shown  on  Fig.  2.6,  2.7  for  various  spatial  locations  z.  At  Fig.  2.6  the 
spectrum  of  X-ray  emission  I(v)  for  various  values  z  is  presented.  One  can  see,  that  the 
spectrum  has  very  stable  form,  practically  not  depending  on  position  z.  It  always  has  the 
distinct  maximum  at  the  energies  Sv  «50  -  60  keV  and  rapidly  falls  down  both  to  the 
smaller  (due  to  photoionization)  and  higher  (due  to  Compton  losses)  values  of  Sy.  Note  that 
this  result  practically  does  not  depend  on  the  value  of  maximal  electric  field  Em.  From  the 
Fig.  2.7  it  is  seen  that  the  maximal  intensity  Imax  of  emission  is  growing  exponentially  with 
electric  field  Em.  The  position  of  emission  maximum  Zmax  is  moving  slowly  in  the  direction 
of  electron  motion. 

The  obtained  results  are  in  a  very  good  agreement  with  all  observations  of  X-ray  spectrum 
during  thunderstorm  at  airplanes  [4],  balloons  [5,6]  and  in  the  mountains  experiments  [6]. 
The  observations  always  show  a  clear  maximum  at  the  energies  ~  60KeV  and  rapid  fall  of 
the  spectrum  both  to  30  KeV  and  100  KeV.  One  can  conclude  from  the  comparison  with  the 
present  theory  that  in  all  these  experiments  the  runaway  avalanche  effect  inside 
thundercloud  was  observed. 

We  note  also,  that  the  effective  growth  of  the  number  of  runaway  electrons  at  EJEc~  1-5-1. 2 
demonstrated  at  Fig.  2.3  lead  to  a  much  stronger  growth  of  the  number  of  low  energy 
electrons  due  to  ionization  of  air  molecules.  These  electrons  rapidly  loose  through  collisions 
with  molecules  their  energy  and  become  thermal  one  [15].  Due  to  attachment  of  the  thermal 
electrons  to  molecules  O2  and  H20  in  a  100  ns  time  scale  the  plasma  of  positive  and  negative 
ions  is  formed.  The  ion  density  in  this  plasma  layer  is  growing  very  effectively  due  to 
runaway  breakdown  process  and  in  a  time  scale  of  the  order  10  s  a  highly  conducting  layer 
could  be  generated.  The  creation  of  this  layer  can  affect  significantly  the  distribution  of  the 
electric  charge  and  lead  to  destabilization  of  thunderstorm  cloud  which  manifests  itself 
eneration  of  lightning.  The  process  of  generation  highly  conducting  region  due  to  runaway 
breakdown  was  discussed  in  details  in  the  previous  papers  [1,7,8]. 


2.5  Conclusion 

Presented  here  theoretical  study  of  runaway  breakdown  process  in  inhomogeneous 
thundercloud  electric  field  has  revealed  two  significant  features  of  X-ray  emission 
generated  by  runaway  electrons: 


1.  The  emission  spectrum  has  well  defined  maximum  located  at  50  -i-  60KeV.  The 
spectrum  falls  down  rapidly  both  to  the  low  energies  (8v  ~  30)KeV  and  to  the  higher 
energies  (Sv  ~  100  -;-  150KeV).  This  statement  practically  does  not  depend  on  z  position, 
where  the  spectrum  is  observed. 

2.  The  intensity  of  X-  ray  emission  is  diminished  effectively  with  the  distance  from  Zmax  - 
position  of  electric  field  maximum  Em=E(zmax),  where  runaway  breakdown  conditions  are 
fulfilled  Em  p  Ec(zmax).  At  the  distance  about  1  km  from  Zmax  the  maximal  intensity  of 
emission  falls  down  more  than  an  order  of  magnitude.  Taking  into  account,  that  according 
to  multiple  observations  of  the  thunderstorm  electric  field  Zmax  is  located  at  the  heights 
4kmp  Zmax  P  8  km  [2,3],  we  deduce,  that  the  X-ray  emission  generated  by  runaway 
avalanche  could  be  observed  at  the  heights  z  p  3km  only. 
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Both  these  statements  are  in  agreement  with  the  existing  observations  [4]-[6],  [16,17].  They 
can  serve  to  a  definite  indication  for  searching  of  runaway  breakdown  effect  in  the 
thunderstorm  atmosphere. 
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Fig.2,1.  The  distribution  function  calculated  for  different  spatial  positions  z.  z=-2  - 
dashed-dotted  line,  z=-l  -  dashed,  z=0  -  dotted  and  z=2.5  -  solid  lines.  Parameter 
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Fig.2.2.  The  electron  distribution  function  F(y,Zmax)  calculated  for  different  values  Em/Ec. 


19 


Fig.2.5.  The  dependence  of  local  X-ray  intensity  lo  generated  by  runaway  electrons  on 
photon  energy  Sy  /o  -  is  a  normalized  flux  proportional  to  jv 
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Fig.2.6.  The  spectrum  of  X  -  rays  calculated  for  different  position  z.  Parameter  £'m/£’c=1.05 
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Fig.2.7.  The  dependence  of  maximal  value  of  X-ray  emission  Imax  and  position  of  this 
maximum  Zmax  on  parameter  Em/Ec. 


3.  Investigation  of  runaway  breakdown  phenomena  in  laboratory  installation  in 

conditions  of  cyclotron  resonance 


Abstract 

The  microwave  electron  cyclotron  resonance  discgarge  in  air  is  sstudied  in  a  magnetic 
mirror  trap.  The  obtained  experimental  data  demonstrate  a  rapid  growth  of  the  number  of 
fast  electrons  with  a  wide  energy  spectrum  (1-300  keV).  The  discharge  begins  with  the 
strong  X-ray  emission,  followed  by  the  syncrotron  emission  at  second  and  third  cyclotron 
harmonics.  The  optic  emission  from  a  low  energy  energy  part  of  electron  spectrum  (<100 
eV)  has  a  delay  up  to  80  ms  with  respect  to  the  X-ray  emission. 

3. 1  Setup  description 

The  setup  represented  in  Fig.3.1  includes  a  vacuum  chamber  2  of  quasispherical 
form,  which  has  a  the  distance  between  a  pair  of  opposite  flanges  450  mm.  The  chamber  is 
placed  between  two  similar  coaxial  copper  coils  1  of  a  static  magnetic  field.  The  magnetic 
field  force  line  distribution  inside  the  chamber  has  configuration  of  a  bottle  with  two 
symmetric  necks  which  forms  a  magnetic  trap  for  collisionless  electrons.  The  trapped 
electrons  confinement  time  in  the  chamber  is  likely  to  be  enough  to  confine  a  bulk  of  initial 
electrons  of  large  energies  during  the  time  of  run-away  electrons  breakdown  in  air.  The 
magnitude  of  the  static  magnetic  field  at  the  center  of  chamber  varies  relative  to  value  875 
Gs  which  corresponds  to  the  ECR  for  the  electric  field  frequency  2,45  GHz. 

A  magnetron  3  generates  electromagnetic  oscillations  at  the  fixed  frequency  2,45 
GHz,  with  varying  power  up  to  5  kW  which  creates  the  ECR  accelerating  field  for  electrons 
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caught  inside  the  mirror  trap.  The  microwave  power  is  transmitted  into  the  chamber  through 
a  waveguide  system  6  which  consists  of  the  adapter  with  pyrex  window  insulating  vacuum 
of  the  chamber  from  athmosphere,  the  horn  anteima  7,  the  circular  waveguide  with  travelling 
TEii  mode,  the  smooth  waveguide  transition  from  a  rectangular  cross  section  to  circular  one, 
the  rectangular  waveguide  elbow,  the  water  attenuator  or  ferrite  circulator  4.  For  purpose  of 
pulse  periodic  regime  realization  it  was  made  mechanical  waveguide  cutter  off.  The 
maximum  microwave  electric  field  strength  at  the  center  of  vacuum  chamber  under 
condition  of  travelling  wave  at  maximum  power  is  EqwISO  V/cm.  In  the  vicinity  of  the  ECR 
the  microwave  electric  field  is  equivalent  to  a  static  electric  accelerating  field  for  electrons. 
Following  experimental  techniques  for  study  of  x-ray  bremmsstrahlung  radiation  and 
optical  radiation  produced  by  the  run  away  electrons  and  microwave  synchrotron  radiation 
at  the  harmonics  of  the  accelerating  field  were  used: 

•  Scintillation  x-ray  transducer  14  (BYCRON,  model  2.5M2.5/3LP-X). 

The  transducer  consists  of  NaJ(Tl)  scintillator  covered  with  aluminium  cup  of  1  mm 
thickness  and  combined  with  a  photomultiplier.  The  optical  axis  of  the  transducer  placed  in 
the  chamber  was  perpendicularly  directed  to  the  trap  axis.  The  operating  quantum  energy 
range  of  the  x-  ray  transducer  100  keV  -5-  1  MeV.  Its  signal  was  registered  at  a  monitor  of 
the  PC  «Pentium»  by  means  of  an  oscillograph  card. 

•  Photomultiplier  16  for  observing  of  visible  air  glow  discharge  under  action  of 
both  ECR-discharge  and  run  away  electrons. 

•  Personal  x-ray  doze  meter  (not  shown)  for  measurement  of  x-ray  spectrum 
outside  the  chamber.  The  operating  quantum  energy  range  of  the  counter  0,05  -s- 
1,5  MeV.  Its  error  of  doze  power  measurements  is  ±  30%;  exposition  time  is  36  s; 
the  dependence  of  doze  measurements  over  the  operating  energy  range  is  no  more 
45%. 

•  Microwave  spectrum  analyser  13  over  the  frequency  range  of  2-nd  and  3-d 
harmonics  for  the  base  frequency  of  the  microwave  generator.  Signals  for  the 
analyser  were  received  from  loop  antenna  12  at  end  of  a  coaxial  feeder 
introduced  along  the  trap  axis. 

•  Detector  1 1  for  monitoring  incident  MW  power. 
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Fig.3.1.  Scheme  of  the  experimental  setup:  1  -  magnetic  system;  2  -  vacuum  chamber;  3  - 
microwave  generator;  4  -  attenuator  or  isolator;  5  -  obturator;  6  -  waveguide 
components;  7  -  horn  coupler;  8  -  vacuum  pump;  9  -  water  attenuator  section;  10  - 
hom  antenna;  1 1  -  detector;  12  -  loop  antenna;  13  -  microwave  receiver  with 
spectrum  analyzer;  14  -  X-ray  receiver;  15  -  optic  fiber;  16  -  photomultiplier;  17  - 
data  display  system  (computers). 

The  insertion  demonstrates  the  magnetic  field  distribution  in  stationary  state 
normalized  to  magnetic  field  in  the  center  of  vacuum  chamber. 
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3.2  Experimental  Results 

In  the  experiments  with  the  described  setup  the  intense  X-ray  emission  was  found  to  be 
generated  under  ECR  discharge  at  the  pressure  range  <10"^  Torr  due  to  bremmsstrahlung 
radiation  of  electrons  accelerated  in  relatively  weak  electric  field  E=5^20  V/cm  up  to 
subrelativistic  energies  100-5-300  keV.  Simultaneously  the  MW  synchrotron  radiation  at  the 
second  harmonic  was  observed.  Its  spectrum  points  out  the  presence  bulk  of  electrons  with 
energies  <10  keV  in  the  system.  Specific  feature  was  that  the  beginning  of  discharge  visible 
emission  owing  to  low  energy  part  formation  of  electron  spectrum  (<100  eV)  had  a  delay 
reaching  80  ms  in  respect  to  the  appearance  of  X-ray  and  synchrotron  signals.  This  allows 
to  separate  the  whole  process  of  the  ECR  breakdown  into  following  two  stages:  first  of 
which  to  be  an  initial  stage  of  breakdown  with  participating  of  high  energy  run  away 
electrons,  and  second  one  to  be  ordinary  ECR  discharge  with  participating  of  low  energy 
electrons.  Such  separation  is  reasonable  due  to  difference  of  physical  processes  resulting  in 
the  characteristic  of  electron  distribution  functions.lt  allows  to  find  analogy  of  ECR 
breakdown  with  the  model  of  runaway  electron  breakdown  by  Gurevich,  Milikh  and 
Roussel-Dupre  applied  for  explanation  of  lightning  mechanism  in  Earth  atmosphere. 
Recently  interesting  results  were  also  observed  in  pulse  mode  of  microwave  power.  Such 
mode  can  be  achieved  by  regular  cutting  off  the  incident  MW  power  with  mechanical 
waveguide  modulator.  The  preliminary  results  are  shown  in  Fig.3.2. 
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One  can  see  (Fig.3.2)  that  at  air  pressure  Torr  we  always  have  a  time  delay 
visible  light  radiation  of  discharge  in  respect  to  the  X-ray  signal  and  the  cut  off  transmitted 
MW  radiation  connected  with  that  process  of  accumulation  of  low  energy  plasma  with 
overcritical  electron  density. 

The  behavior  of  signals  in  Fig.3.3  corresponding  to  air  pressure  2-10'^  demonstrates 
the  visible  discharge  is  ahead  in  respect  to  the  X-ray  signal.  However  one  can  see  the 
flickering  of  discharge  when  the  incident  microwave  energy  is  redistributed  between 
runaway  electrons  and  low  energy  plasma. 

The  main  objectives  of  our  future  work  are  detailed  investigations  of  runaway 
electron  breakdown  taking  into  account  its  temporal  and  spatial  dynamics,  definition  of 
energy  characteristics  of  the  process  and  energy  balance,  measurements  of  energy  y-quanta 
spectra  depending  on  experimental  conditions. 

To  study  the  dynamics  of  the  discharge  plasma  density  due  to  the  RB  a  microwave 
interferometer  (A,=8  mm)  of  opened  resonator  type  will  be  constructed.  The  x-ray 
spectrometer  in  the  energy  range  30  keV  -  1  MeV  will  be  created  with  use  of  the  X-ray 
scintillation  detector  (in  our  possession)  and  the  pACE™-system  (EG&G  ORTEC)  which 
consist  of  amplifier,  multichannel  analyzer  plug-in  card  with  emulation  software(will  be 
ordered  to  the  ISTC).  The  averaged  x-ray  spectra  depending  on  the  regimes  of  the  RB  and 
optical  spectra  of  the  RB  plasma  will  be  studied.  Both  probe  measurements  and  electron 
energy  distribution  function  measurements  with  an  electrostatic  analizer  for  electrons 
leaving  the  magnetic  trap  will  be  carried  out. 

We  intend  to  continue  investigation  with  pulse  mode  of  MW  generation  because  it’s 
reasonable  to  have  the  time  limited  action  on  the  RB  for  study  of  its  dynamics.  For  this 
purpose  the  waveguide  modulator  5  (Fig.3.1)  cutting  off  continuous  microwave  radiation 
was  constructed.  It  allows  to  vary  duration  and  period  of  pulses  in  wide  range.  The  demands 
for  the  modulator  was  caused  by  the  reason  that  MW  generator  is  able  to  work  in 
continuous  mode  generation,  but  not  able  to  do  it  in  pulse  mode  if  to  use  pulse  high  voltage 
supply.  Microwave  pulses  after  the  modulator  have  trapezoidal  form  with  rise  and  fall  time 
of  order  1/3  of  plateau.  Owing  to  the  pulse  mode  we  hope  to  study  lifetime  of  high  energy 
electrons  caught  by  the  magnetic  trap  during  a  previous  microwave  pulse  and  their  role  as 
seed  in  the  RB  phenomenon  during  followed  pulse. 


28 


4.  Investigation  of  an  influence  of  combined  effects  of  runaway  breakdown  and 
cosmic  rays  on  lightning  processes  in  thunderstorm  atmosphere 

4. 1  Abstract 

The  intensive  X-ray  emission  during  thunderstorms  was  studied  at  the  height  3340  m  above 
sea  level.  For  the  first  time  the  short-time  (1 — 5  min)  bursts  of  X-ray  emission  were 
observed.  The  bursts  are  highly  correlated  at  the  wide  space  region  (about  0.5  km).  The 
main  component  of  X-ray  emission  in  bursts  has  the  energy  50 — 80  keV.  The  observed 
bursts  could  be  attributed  to  the  bremsstrahlung  determinated  by  the  runaway  breakdown 
(RB)  effect  in  the  thundercloud  electric  field.  Ground  based  observations  of  RB  open  a  wide 
opportunities  for  the  studies  of  fundamental  processes  in  thunderstorms. 

4.2  Introduction 

The  new  mechanism  of  lightning  initiation  in  thunderstorm  clouds  based  on  a 
combined  action  of  runaway  breakdown  (RB)  and  extensive  atmospheric  showers  (EAS) 
was  proposed  in  [1].  According  to  this  theory  in  shower  core  the  ionization  of  the  air  due  to 
RB  effect  is  growing  strong  enough  to  produce  a  local  highly  conductive  plasma.  The 
thunderstorm  electric  field  is  amplified  strongly  inside  such  region.  That  can  serve  to  the 
triggering  of  a  lightning  discharge. 

The  Lebedev  Physics  Institute  (LPI)  group  proposed  an  experimental  program  to 
study  this  mechanism.  The  program  includes  simultaneous  observations  of  cosmic  ray 
showers,  runaway  breakdown  emission  and  lightning  discharges  at  the  mountain  cosmic  ray 
station  (MCRS).  The  LPI  MCRS  is  situated  on  a  pass  in  Tien-Shan  mountains  at  the  height 
3340  m  above  sea  level  (effective  atmosphere  depth  690  g/cm^).  The  thunderstorm  clouds 
go  through  the  pass  at  the  height  about  0.1 — 1  km  during  of  three  months  yearly.  The 
installation  can  detect  extensive  air  showers  measuring  simultaneously  EAS  electrons  with 
energies  0.1-h10  MeV,  muons  with  energies  above  5  GeV  and  neutrons  with  energies  above 
100  MeV.  The  RB  effect  lead  to  a  strong  amplification  of  the  X-ray  emission  inside  and 
below  the  thunderclouds  [2]-[6].  To  observe  RB  emission  a  system  of  Geiger  counters  was 
installed  in  a  wide  region  (approximately  500  m  around  MCRS).  To  fix  the  position  of 
lightning  discharge  by  its  strong  electromagnetic  and  acoustic  radiation  a  group  of  special 
antennas  together  with  directed  microphones  was  established. 

In  the  present  work  we  inform  about  the  first  results  of  X-ray  emission 
measurements.  The  strong  amplification  of  emission  during  the  whole  thunderstorm  period 
is  found.  For  the  first  time  the  bursts  —  short  time  strong  peaks  of  X-ray  radiation  were 
observed.  The  typical  time  length  of  such  a  burst  is  1^5  min  and  radiation  increases  are 
highly  correlated  in  different  points  at  the  scales  minimum  of  the  order  of  several  hundred 
meters. 

4.3  Experimental  Installation 

The  schematic  layout  of  MCRS’s  detectors  disposition  is  shown  at  Fig.4.1.  The 
installation  used  in  this  experiment  consists  of  the  following  main  parts: 

•  The  shower  system  on  the  base  of  scintillation  detectors  located  nearly  concentric 
circumferences  on  the  flat  part  of  the  pass  (the  “Hadron”  experiment  installation  [7]). 
The  scintillators  can  detect  the  high  energy  (above  10  MeV)  electrons  and  gamma  quanta 
travelling  through  the  installation  within  an  extensive  air  shower  fronts. 

•  The  NM64  type  neutron  supermonitor  [8]  consisting  of  three  independant  units.  Each 
unit  has  3  m^  sensitive  area  and  includes  six  registration  channels  on  the  base  of  neutron 
counters  SNM15.  The  supermonitor  was  designed  for  the  registration  of  one-minute 
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intensity  variations  of  cosmic  ray  hadrons  with  energy  above  50 — ^100  MeV.  The 
electronics  of  neutron  monitor’s  registration  channels  is  designed  to  provide  methodic 
error  at  the  level  of  0.01%  in  intensity  measurements.  The  input  threshold  value  for 
neutron  counter's  pulses  is  set  to  be  equal  to  10  mV,  corresponding  background  pulse 
intensity  being  2700  pulses  per  min  (p.p.m.)  for  each  unit. 

•  The  boxies  of  high  threshould  gamma  ray  and  electron  detector  are  placed  above 
supermonitor  units  as  shown  at  Fig.4.  lb.  This  detector  is  made  on  the  base  of  Geiger- 
Mueller  counter  tubes  of  SI5G  type  which  are  capable  to  register  hard  X-  and  soft  y-ray 
emissions  with  nearly  1%  efficiensy.  Four  detector  boxies  each  of  wich  contains  20 
counters  with  parallel  outputs  and  corresponds  to  the  one  registration  channel  were 
mounted  above  each  of  three  supermonitir  units,  their  total  sensetive  area  being  equal  to 
8  m^.  Because  the  counters  of  this  detector  are  operating  in  proportional  regime  it  is 
necessary  to  decrease  their  dead  times  up  to  5-10  mcs  (which  is  20  times  shorter  as  for 
usual  Geiger  counters).  To  achive  the  charged  particles  registration  efficiency  as  high  as 
for  Geiger  regime  (above  95%)  the  registration  chaimel's  threshould  is  set  to  20  mV. 
The  roof  of  die  building  where  neutron  supermonitor  and  high  threshould  y-ray  detector 
is  situated  is  made  of  0.8  mm  thick  iron,  the  wood  ceiling  has  thikness  about  15  cm  and 
the  detector’s  boxes  are  made  of  2  mm  thick  aluminium.  It  provides  the  upper  threshold 
about  25 — 30  MeV  for  relativistic  charged  particles  registration.  The  main  absorber 
above  the  detector  is  carbon  with  the  effective  thikness  13.6  g/cm^.  Thus,  the  threshold 
value  for  the  photons  is  about  30  keV.  The  general  count  intensity  for  all  channels  is 
determinated  mainly  by  the  electron  background  flux  and  is  equal  to  3.2T0^  p.p.m.  The 
corresponding  statistical  error  in  1 -minute  data  is  less  than  0.2%. 

•  The  muon  detector  made  of  SI5G  type  counters  with  common  sensitive  area  2.7  m^  is 
placed  under  one  of  the  supermonitor  units  (see  Fig.4.  lb).  Output  pulses  from  each  of 
four  lower  detector  boxies  are  switched  to  a  coincidence  scheme  with  the  pulses  coming 
from  corresponding  boxies,  located  above  supermonitor  unit.  The  thickness  of  matter 
above  lower  boxies  is  15  g/cm^  of  carbon  and  150  g/cm^  of  lead  absorber.  The  latter 
exeeds  35  electromagnetic  cascade  units,  excluding  completly  the  registration  of 
atmospheric  photons  and  electrons  with  energies  below  10  GeV.  The  detector  registrates 
mainly  the  muons  with  energies  above  220  MeV  travelling  through  the  unit  in  nearly 
vertical  direction.The  general  count  intensity  is  about  7-10'*  p.p.s.  Hence,  the  fluctuations 
in  one-minute  data  are  about  0.5%.  The  registration  channel’s  characteristics  of  muon 
detector  are  just  the  same  as  those  for  the  high  threshould  gamma  detector. 

•  The  newly  build  low  threshould  wide  spread  X-ray  detector  (WS)  created  in  Summer 
1999  is  designed  to  investigate  the  thunderstorm  clouds  low  energy  electromagnetic 
emission.  It  consists  of  four  registration  points  situated  perpendicular  to  the  preferential 
direction  of  thunderstorm  clouds  motion  over  the  cosmic  ray  station's  pass  (see  Fig.4.  la). 
The  distance  between  the  extreme  right  and  extreme  left  points  is  about  500  m.  Each  of 
the  points  has  two  layers  of  sixty  SI5G  counters.  The  outputs  of  the  counters  in  each 
layer  are  combined  in  one  registration  channel.  Hence,  there  are  two  registration 
channels  in  each  point  with  sensitive  square  2  m^  each.  The  SI5G  counters  are  put  in 
Geiger-Mueller  regime,  and  the  threshold  of  registration  channels  is  equal  to  2 — 3  V. 
Electronic  circuits  in  each  detector  point  are  based  on  vacuum  tubes  in  order  to  protect 
them  from  the  nearby  lightning  discharges.  The  powerfull  vacuum  tubes  ensure  the 
output  pulse  front  length  not  worse  then  2 — 3  mcs  at  output  resistance  of  240  Ohm  and 
capacity  about  20  mcF.  Due  to  this  distortion  of  the  pulse  form  during  its  passage 
through  300 — 400  m  long  cable  is  negligible.  The  detectors  of  each  registration  points 
are  placed  within  a  little  light  house.  The  counters  are  separated  from  the  atmosphere  by 
the  0.8  mm  (or  0.74  g/cm^)  thick  iron  roof 
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All  the  detectors  described  above  have  the  bouth  independed  power  supply  and  data 
recording  systems.  Distant  detectors  are  connected  to  a  central  registration  system  (located 
in  the  building  with  neutron  supermonitor)  by  screened  underground  cables.  The  data  are 
transfered  to  computer  by  the  means  of  a  CAMAC  system.  To  avoid  the  danger  of 
semiconductor  electronics  damage  in  the  moment  of  lightning  discharge,  the  cables  are 
conjugated  with  CAMAC  scaler  blocks  throughout  vacuum  tube  based  circuits. 

The  pulses  coming  from  all  of  eight  channels  of  low  threshould  WS  X-ray  detector 
were  counted  by  CAMAC  scalers  during  of  one  second  time  intervals  and  accumulated 
numbers  of  them  were  recorded  by  registration  computer.  Amounts  of  pulses  obtained  from 
other  detectors  were  counted  during  of  one  minute  intervals  and  registrated  together  with 
the  clock  indication,  temperature  and  barometric  pressure.  The  registration  was  perfomed 
continuously  during  of  4  monthes  in  July — October  1999. 

We  have  tested  the  electric  noise  in  all  channels  by  switching  off  the  high  voltage 
power  supply  in  registration  points.  In  such  a  case  the  pulses  from  the  channels  were  not 
registrated  even  in  thunderstom  conditions. 

4.4  Experimental  results 

Intensity  variations.  During  of  July  and  up  to  the  last  days  of  August,  1999  there 
were  no  thunderstorms  in  the  vicinity  of  mountain  station.  The  typical  measurements  of 
pulse  intensity,  obtained  from  various  detectors  at  this  quiet  period  are  presented  at  the  3.2. 
The  shown  data  correspond  to  the  total  intensity  registrated  by  all  channels  of  each  detector 
(i.e.  the  amounts  of  pulses  obtained  from  8  channels  of  low  threshould  X-ray  detector,  12 
channels  of  high  threshould  y-ray  detector,  3  neutron  monitor  units,  4  muon  coincidence 
channels  show  the  intensity  curves  at  this  figure). 

One  can  see  that  both  low  and  high  threshould  y-ray  detectors  exhibit  no  significant 
changes  in  count  velocity  through  the  whole  period  shown  (6  houres).  The  signals  from 
neutron  monitor  and  muon  detector  demonstrate  the  same  stability  (their  slow  trend  is 
caused  by  variation  of  atmospheric  pressure). 

In  the  end  of  August  1999  we  have  observed  five  cases  when  at  low  heights  over  the 
MCRS  were  passing  the  thunderstorm  cloud  fronts  accompanied  by  the  strong  lightning 
discharges.  The  same  data  as  given  at  Fig.4.2  for  a  quiet  day  are  presented  at  Fig.4.3a-e  for 
typical  thunderstorms  events.  Time  scale  at  the  bottom  of  the  figure  is  centered  relatively  to 
the  moment  when  thunderstorm  beginning  was  fixed.  This  moment  is  marked  by  a  vertical 
dotted  line.  Below  we  shall  discuss  the  typical  thunderstorm  event  on  August,  31. 
Corresponding  data  are  presented  at  Fig.  4.3.e. 

It  can  be  seen  from  the  typical  curves  shown  at  Fig.4.3e,  that  each  time  after  the 
moment  of  thunderstorm  beginning  a  dramatic  growths  of  emission  intensity  measured  both 
by  low  and  high  threshould  y-ray  detectors  occures.  We  emphasize  the  significant  difference 
between  these  two  types  of  detectors:  low  threshold  detectors  measured  practically  free  flux 
of  X-rays,  and  high  ftneshold  detector  that  was  under  the  roo/having  additional  absorption 
equivalent  to  13.6  g/cm^  of  carbon.  One  can  see  from  the  Fig.4.3.e,  that  flie  signals  of  the 
free  (low  threshould)  and  under  the  roof  (high  threshould)  detectors  are  well  correlated  in 
time.  These  correlations  are  quite  regular  having  the  same  character  for  every  thunderstorm. 
The  characteristic  scale  of  long  time  emission  increses  is  2 — 3  houres  after  the  moment 
when  thunderstorm  clouds  begin  to  pass  the  MCRS. 

At  the  same  time  the  neutron  monitor  and  muon  detector  (which  are  sensitive  to 
energetic  particles  of  cosmic  radiation  only)  do  not  demonstrate  any  thunderstorm  caused 
effect. 

The  results  obtained  for  the  all  of  five  observed  thunderstorms  are  quite  analogous. 
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Space  correlations.  The  y-ray  intensity  increases  during  thunderstroms  discussed 
above  have  been  studied  in  data  streams  obtained  separately  from  each  of  low  threshould 
detector  points  shown  at  Fig.4.1a.  The  correlation  coefficients  between  all  available  pairs  of 
detector  points  were  determined.  The  results  —  time  dependences  of  correlation  coefficients 
Cij  calculated  for  100  sec.  Length  time  periods  for  three  detector  points  pairs  —  are 
presented  at  Fig.4.4  for  one  of  the  observed  thunderstorms  (the  corresponding  data  obtained 
for  all  other  detector  pairs  and  of  other  thunderstorms  are  quite  analogous).  One  can  see 
from  the  Fig.4.4  that  under  the  fine  weather  conditions  the  data  at  different  observational 
points  are  practically  noncorrelated.  During  the  thunderstorm  period  the  picture  is 
absolutely  different:  correlation  coefficients  for  all  low  threshold  detector  points  pairs 
becomes  very  high  —  about  0.8  and  even  higher.  This  high  correlation  demonstrate  that 
the  observed  of  X-ray  intensity  strong  growth  is  deeply  connected  with  thunderstorm 
processes  developing  simultaneously  at  the  wide  space  region  of  the  order  of  0.5  km,  maybe 
even  larger. 

Short  time  emission  bursts.  The  accompanying  thunderstorm  emission  intensity 
variations  are  well  correlated  in  time  also.  It  can  be  seen  from  Fig.4.3.e  for  X-ray  signals  of 
the  free  (curve  1)  and  under  the  roof  (curve  2)  detectors.  Analogous  time  eorrelations  exist 
for  the  emission  at  different  low  threshold  detector  points.  These  eorrelations  are  quite 
regular,  they  have  the  same  character  for  every  thunderstorm.  The  characteristie  timescale 
of  this  longtime  emission  variation  is  30-^60  min.  One  can  mention  the  following  peculiarity 
of  emission  intensity  —  it  has  often  a  double  or  triple  pike  structure  (see  Fig.4.3.e). 

An  interesting  new  feature  is  short  time  bursts  of  X-ray  emission.  Characteristic 
length  of  the  bursts  is  1-5-5  min.  The  typical  short  time  bursts  occured  during  of 
thunderstorm  event  of  31  August  1999  can  be  seen  at  Fig.4.3.e  in  time  moments  0.2  (peak 
1),  0.4  (peak  2),  1.4  (peak  3)  houres  after  the  beginning  of  thunderstorm.  It  is  obvious  good 
correlation  between  the  intensity  bursts  of  low  and  high  threshould  emission  detectors.  Any 
increases  at  corresponding  time  moments  in  the  intensity  of  energetic  cosmic  ray 
components  (i.e.  neutrons  and  muons)  are  absent. 

The  one  of  the  bursts  data  for  low  threshould  detector  intensities  (peak  3)  are  shown 
in  large  scale  at  Fig.4.5  both  for  the  sum  of  all  channels  and  for  each  detector  point 
separately.  One  can  conclude  from  this  figure,  that  emission  intensity  bursts  are  correlated 
strongly  between  various  detectors  separated  with  space  distancies  of  the  order  of  some 
hundreds  of  meters.  According  to  statistical  errors  marked  at  Fig.4.5,  the  value  of  intensity 
increase  during  of  a  typical  emission  burst  exeeds  the  level  of  statistical  flluctuations  by 
15 — 30  times. 

Absorption  of  emission.  The  two  types  of  measurements  of  emission  intensity  were 
performed  simultaneously:  low  threshold  and  high  threshold,  which  demonstrate  quite  a 
good  time  correlation  (Fig.4.3.e).  The  following  main  difference  between  these 
measurements  exists.  The  low  threshold  detectors  were  separated  from  the  atmosphere  by 
the  thin  (0.8  mm)  iron  protector  only.  The  high  threshold  detectors  had  the  same  protector, 
but  they  were  installed  inside  flie  neutron  supermonitor  building,  which  has  a  roof  with 
effective  thickness  13.6  g/cm^  of  carbon.  The  roof  gave  a  significant  additional  absorption. 
By  the  comparison  of  the  data  from  free  (low  threshold)  and  under  the  roof  (high  threshold) 
detectors  the  characteristic  absorption  of  emission  could  be  found. 

To  perform  comparison  we  determined  first  five-minute  averaged  background  of 
both  sets  of  data.  As  fast  electrons  determine  the  major  part  of  the  counts  we  have  to 
eliminate  the  background  and  compare  the  excess,  which  is  due  to  X-ray  effect.  The  excess 
relation  r  for  both  sets  of  data  is  presented  at  Fig.4.6.  One  can  see,  that  the  typical  value  of 
excess  relation  of  the  free  (low  threshold)  to  the  under  roof  (high  threshold)  data  is  2-5-4. 
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Such  an  absorption  in  13.6  g/cm^  carbon  is  characteristic  for  energy  of  y-quanta  200-5-500 
keV[18]. 

It  is  important,  that  during  short  time  bursts  the  value  of  r  is  always  significantly 
higher,  growing  up  to  8-5-12  (Fig.4.6).  That  shows  that  the  energy  spectrum  of  X-rays  during 
short  bursts  is  essentially  different  from  the  main  part  of  thunderstorm:  its  major  component 
has  low  energy  X-rays  (<  50-5-80keV)  which  are  effectively  absorbed  by  the  roof  Note  that 
peaks  are  more  pronounced  in  absorption  curves  (Fig.4.6),  than  in  intensity  (Fig.4.3.e).  It 
opens  a  possibility  to  use  absorption  as  an  effective  method  of  X-ray  burst  detection. 

Flux  of  X-ray  quanta.  To  determine  the  flux  of  X-ray  quanta  we  have  first  to  average 
the  excess  counts  of  low  threshold  data  and  then  takes  into  account  absorption  by  the  iron 
protector  and  sensitivity  of  Geiger-Muller  counters  SI5G  for  X-ray  quanta.  That  gave  an 
approximate  value  of  flux  in  bursts  10^20  photons/(cm  sec  sr). 

4.5  Discussion 

The  aim  of  this  work  is  to  study  the  influence  of  thunderstorms  on  the  X-ray  emission  from 
atmosphere.  Effective  growth  of  the  emission  in  all  observed  thunderstorms  is  definitely 
established. 

In  fine  weather  conditions  weak  disturbances  of  the  emission  (~  1%)  determined  partly  by 
the  changes  in  atmospheric  pressure  are  regularly  observed  (Fig.4.1)  That  agrees  well  with 
the  previous  statements  [9]  -  [11].  During  thunderstorm  the  picture  changes  fundamentally. 
The  X-ray  emission  is  strongly  amplified  for  hours.  The  emission  amplification  in  the 
thunderstorm  period  was  studied  previously  in  a  number  of  papers  (see  [10]-[12]). 

The  main  peculiarities  of  our  observations  are: 

1.  X-ray  emission  was  studied  at  the  ground,  but  high  in  the  mountains  (height  3340m). 
Analogous  investigations  were  done  previously  only  by  Shaw  [10]  and  Suszcynsky  et 
al.  [1 1]  (mainly  at  lower  heights). 

2.  For  the  first  time  the  emission  was  observed  simultaneously  in  several  points  with  the 
large  distance  between  them  (Fig.4.1). 

3.  For  the  first  time  the  free  and  partly  screened  emission  was  studied  simultaneously. 

Below  we'll  compare  our  observations  with  the  previous  studies.  The  new  results  will  be 
discussed. 

The  main  mechanisms  of  strong  amplification  of  X-ray  emission  intensity  during 
thunderstorm  are: 

1 .  Gamma  ray  emission  from  the  beta  decay  of  radon  daughter-ion  products. 

2.  Bremsstrahlung  of  accelerated  electrons  due  to  the  action  of  a  strong  thunderstorm 
electric  field. 

Radon  effect  is  determined  by  the  precipitation  of  radon  and  its  products  to  the  ground 
during  rainfall.  Gamma  ray  emission  from  radon  daughter-ion  P  decay  has  energies  of 
quanta  0.2-5-2  MeV.  The  decay  processes  and  rainfall  determine  its  main  time 
characteristics:  20  -5-  30  min  growing  up  and  about  50  min  after  washout  [11,13]. 

Electric  field  effect  is  determined  by  the  acceleration  of  electrons.  Wilson  [14]  was  the  first 
who  suggested  that  the  secondary  electrons  generated  by  cosmic  rays  could  be  accelerated 
by  the  thunderstorm  electric  field.  A  new  fundamental  step  was  done  in  the  runaway 
breakdown  (RB)  theory  developed  by  Gurevich  et  al.  [2].  It  was  shown  that  accelerated 
electrons  could  organize  an  exponentially  growing  avalanche  in  thunderstorm  electric  field. 
RB  effect  is  supposed  to  determine  maximal  electric  field  reached  in  thunderclouds 
[2,15,16].  Directly  RB  process  we  discussed  in  the  Introduction  in  connection  with  the 
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lightning  triggering  theory  [1].  RB  effect  lead  to  flie  strong  amplification  of  X-ray  emission 
due  to  the  bremsstrahlung. 

Amplification  of  X-ray  emission  attributed  to  RB  effect  was  observed  in  [3,4,17]  during  the 
measurements  inside  thunderclouds  at  the  balloons.  The  main  features  of  these  observations 
could  be  summarized  as  follows: 

a)  The  intensity  of  X-ray  emission  at  the  height  ~  4  km  is  sharply  (in  seconds  time  scale) 
growing  up  to  10^  in  the  comparison  with  the  standard  level  determined  by  cosmic  rays. 
The  maximum  of  observed  emission  bursts  lays  in  the  diapason  60-5-90  keV,  the 
maximal  count  rate  is  40-5-50  photon/(cm^  sec  sr)  [4,17]. 

b)  The  horizontal  scale  of  the  region  of  emission  amplification  is  of  the  order  1-km  or 
more  in  agreement  with  the  layered  structure  of  thunderstorm  clouds. 

c)  The  time  scale  of  X-ray  bursts  is  of  the  order  1  min. 

d)  The  intensity  of  bursts  diminishes  sharply  (in  second  time  scale)  with  the  lightning.  But 
in  a  10-sec  time  scale  it  can  grow  up  again. 

e)  Lightning  was  seen  during  all  bursts.  So,  it  is  natural  to  suppose,  that  a  correlation  of 
lightning  events  with  X-ray  intensity  bursts  do  exist. 

Let  us  discuss  now  results  of  our  observations,  taking  into  account  the  main  features  of 
both  considered  mechanisms  of  X-ray  emission  amplification. 

1.  The  general  structure  of  a  long  time  emission  (Fig.4.3.e)  agrees  well  with  the  radon 
source.  A  peculiarity,  which  could  be  mentioned,  is  the  usual  double  pike  structure  of 
our  intensity  curves  (Fig.4.3.e).  The  energy  of  gamma  quanta  (0.2-5-0.5  MeV), 
determined  by  the  absorption  in  the  13.6  g/cm^  carbon  for  the  main  part  of  observed 
thunderstorm  emission  (Fig.4.6),  agrees  with  radon  source  also. 

2.  Short  time  bursts  of  X-ray  emission  (Fig.4.3.e-4.6)  are  the  main  new  result  obtained  in 
our  measurements.  The  bursts  are  well  correlated  at  the  scales  0.5-km  or  more.  Such 
correlated  short  time  bursts  in  the  ground-based  measurements  are  observed  for  the  first 
time.  There  exists  a  quite  good  correlation  between  the  observations  of  free  and  partially 
absorbed  bursts.  The  analysis  of  absorption  effect  shows,  that  the  energy  spectrum  of  X- 
rays  in  the  bursts  is  essentially  different  from  the  main  part  of  thunderstorm:  its  major 
part  consists  of  the  low  energy  X-rays  (E  <  50-5-80  keV).  Such  type  of  energy  spectrum 
is  characteristic  for  RB  effect  only. 

In  Table  1  the  main  features  of  our  observations  of  bursts  are  compared  with  the  X-ray 
emission  measured  at  balloons.  We  see  that  the  features  of  the  balloon  measurements 
attributed  to  RB  effect  are  rather  close  to  the  short  burst  observed  on  the  ground  in  our 
measurements  (Table  1).  One  can  conclude  that  in  short  bursts  the  runaway  breakdown 
bremsstrahlung  was  observed. 

We  note  that  this  statement  does  not  contradict  to  the  other  ground  based  studies,  including 
the  mostly  comprehensive  and  detailed  one  made  by  Suschinskhy  et  al.  [1 1].  As  in  [1 1],  the 
main  part  of  emission  in  our  study  definitely  belongs  to  the  radon.  In  the  possibility  to 
observe  RB  emission  effect  in  our  case  essential  role  could  be  played  by  higher  sensitivity 
of  our  system  and  observations  of  X-ray  absorption.  One  should  mention  also  a  significant 
difference  in  heights:  the  main  part  of  measurements  in  [11]  were  done  at  the  height  2000 
m,  and  our  —  at  3340  m.  We  emphasize  that  the  mostly  probable  region  of  the  high  values 
of  electric  field  which  can  lead  to  RB  effect  in  thunderstorm  clouds  E~Ec(z)  according  to 
multiple  balloon  measurements  [15,16]  lays  at  the  heights  4-5-8  km.  The  characteristic  length 
of  energetic  photon  E  ~  30  -5-  100  keV  propagation  in  atmosphere  is  If  <1  km  [6].  So,  the 
probability  to  observe  RB  emission  is  strongly  diminishing  with  height  below  3.5  km. 

We  see  that  though  the  RB  emission  is  effectively  screened  by  radon  it  could  be  brought  to 
light  in  the  ground  based  measurements  as  well.  Of  course  more  detailed  ground  based 
studies  including  higher  heights  and  better  energy  spectrum  and  time  resolution  are  needed. 
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4.6  Conclusion 

In  conclusion  we  will  formulate  briefly  the  main  results  of  our  observations. 

1 .  X-ray  emission  during  thunderstorm  in  mountains  at  the  height  3340  m  is  studied  for  the 
first  time  simultaneously  in  several  different  places. 

2.  It  is  found  that  the  emission  is  well  correlated  at  the  scale  of  the  order  0.5  km  or  more. 

3.  For  the  first  time  well  defined  short  bursts  of  X-ray  emission  are  determined.  The  bursts 
are  highly  correlated  at  the  whole  studied  space  region  (~0.5km). 

4.  The  noticeable  diminishing  of  emission  intensity  penetreating  through  the  13.6  g/cm^ 
carbon  layer  is  established.  It  means  that  the  main  component  of  X-rays  observed  in 
bursts  has  low  energy  (50-5-80  keV). 

5.  The  X-ray  flux  observed  in  bursts  could  be  estimated  as  10-5-20  photon/(cm^  sec  sr). 

The  preliminary  analysis  of  our  measurements  allows  to  suppose  that  the  observed 
bursts  of  X-ray  emission  is  the  result  of  bremsstrahlung  amplification  determined  by  the 
runaway  breakdown  effect  in  thundercloud  electric  field. 
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Table  1 


Ground 

Balloon 

Time  duration 

1  ^  5  min 

>  1  min 

Space  correlation 

>  0.5  km 

>  1  km 

Characteristic  energy 

50-80keV 

60  ^  90  keV 

Characteristic  intensity 

10  20  photons/cm^  sec  sr 

40  -i-  50  photons/cm^  sec  sr 

Heights 

3340  m 

3500  4000  m 
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GAMMA  RAY  DETECTOR  BOX: 

20  SI5G  TYPE  COUNTERS 
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Figure  4.1 .  Disposition  layout  of  mountain  cosmic  ray  station’s  detectors  (A):  1-  the 
center  of  shower  scincillation  system,  2-  neutron  monitor  and  high-threshould 
gamma  detectors,  3  -  the  system  of  four  WS  low-threshould  gamma-  and  X-ray 
detector  points,  (B):  the  cross-section  of  a  neutron  monitor  unit  together  with  high- 
threshould  and  muon  detectors. 
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Figure  4.2.  Typical  behaviour  of  various  cosmic  ray  components  and  atmospheric 
parameters  in  the  absence  of  thunderstorm. 
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Figure  4.3a.  The  behaviour  of  various  cosmic  ray  components  and  atmospheric 
parameters  during  of  the  thunderstorm  at  26  August  1999. 
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Figure  4.3b.  The  behaviour  of  various  cosmic  ray  components  and  atmospheric 
parameters  during  of  the  thunderstorm  at  27  August  1999. 
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Figure  4.3c.  The  behaviour  of  various  cosmic  ray  components  and  atmospheric 
parameters  during  of  the  thunderstorm  at  29  August  1999. 
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Figure  4.3cl.  The  behaviour  of  various  cosmic  ray  components  and  atmospheric 
parameters  during  of  the  thunderstorm  at  29  August  1999. 
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Figure  4.3.e.  The  behaviour  of  various  cosmic  ray  components  and  atmospheric 
parameters  during  the  thunderstorm  on  31  August  1999. 
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Figure  4.4.  Time  dependences  of  correlation  coefficients  Cy  for  three  detector  pairs 
during  the  thunderstorm  on  31  August  1999. 
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where  each  of  i  and  j  indices  designates  one  of  the  low  threshould  detector  points 
(if  j=1^4,  i#),  ni  and  nj  mean  amounts  of  pulses  obtained  from  bouth  registration 
channels  of  i-th  and  j-th  detector  during  of  1  sec  intensity  counting,  <ni>  and  <nj> 
are  the  corresponding  average  values  defined  over  the  measurements  immediately 
preceding  thunderstorm  moment.  Summation  in  the  above  formulae  was  carried 
out  during  of  100  second  length  time  periods. 


44 


Figure  4.6.  The  excess  relation  of  low  threshold  to  high  threshold  data 
during  the  thunderstorm  on  31  August  1999. 
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Studies  of  Intensive  X,  y  and  Optic  Emission  in  Runaway 

Breakdown  Process 

This  work  is  driven  by  intersection  of  two  lines  of  activities  in  studying  of  the 

electromagnetic  emission  generated  by  the  suprathermal  electrons.  We  have  addressed  the 

following  advanced  problems  in  the  grant  proposal: 

1 .  Theoretical  investigation  of  runaway  breakdown  process  in  stationary  electric  field. 

2.  Investigation  of  runaway  breakdown  phenomena  in  laboratory  installation  in  conditions 
of  cyclotron  resonance. 

3.  Investigation  of  an  influence  of  combined  effects  of  runaway  breakdown  and  cosmic 
rays  on  lightning  processes  in  thunderstorm  atmosphere. 
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Shepetov  A.L.,  Vildanova  L.I.,  Zybin  K.P.,  Gurevich  A.V.  "Intensive  X-ray  Emission 
Bursts  During  Thunderstorm" 
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1.  Theoretical  investigation  of  runaway  breakdown  process  in  stationary  electric 

fieid 


1. 1  Abstract 

The  full  kinetic  theory  of  runaway  breakdown  for  electrons  at  high  energy  range  is 
developed.  The  theory  takes  into  account  both  ionization  and  bremsstrahlung  processes.  The 
electron-positron  pairs  generation  is  studied.  The  possibility  to  observe  e^e”  annihilation 
line  during  intensive  discharges  in  atmosphere  is  predicted. 

12  Introduction 

Runaway  breakdown  is  a  new  phenomena  determined  by  the  effect  of  constant 
electric  field  at  high  energy  (c  >  100  Kev)  electrons  [1].  Its  physical  mechanism  is  deeply 
related  to  the  Coulomb  character  of  neutral  molecules  ionization  by  high  energy  electrons: 
collisional  crossection  a  in  nonrelativistic  region  drops  rapidly  with  energy  increase  a  oc 
lls^ .  Thus  fast  electrons  in  the  presence  of  electric  field  E>  Ec  {Ec  -  critical  field)  are  able 
to  become  runaways:  in  other  words  they  can  be  accelerated  by  the  electric  field. 
Accelerated  electrons  produce  a  number  of  low  energy  electrons  by  ionization  of  neutral 
molecules.  In  addition,  a  small  number  of  high  energy  electrons  may  be  generated.  These 
electrons  become  runaways  also  if  their  energy  exceeds  critical  value  Sc'^c^(Ec/E)  .  Being 
accelerated,  they  produce  another  generation  of  fast  runaway  electrons  and  the  process  is 
ongoing,  i.e.  as  a  result  an  avalanche  of  runaway  electrons  starts.  Exactly  this  process  is 
called  "runaway  breakdown".  It  is  determined  by  multiplication  and  acceleration  of  fast 
electrons.  The  exponential  growth  of  the  runaway  electrons  number  characterizes  the 
development  of  the  avalanche  and  absorption  of  electromagnetic  energy  by  discharge. 

The  kinetic  theory  of  runaway  breakdown  was  developed  previously  [2-5].  In  these 
papers  electrons  with  not  very  high  energies  (f  <  5  MeV)  were  considered.  The  goal  of  the 
present  work  to  provide  the  full  kinetic  description  of  runaway  breakdown  effect  for  high 
energies  also.  We  will  demonstrate  that  at  high  energies  the  main  role  in  runaway 
breakdown  process  is  played  not  by  ionization  Coulomb  collisions,  but  by  much  more 
stronger  bremsstrahlung  effect.  The  structure  of  electron  distribution  function  at  high 
energies  would  be  determined.  The  role  of  energy  distribution  of  initial  fast  electrons 
triggering  the  runaway  avalanche  will  be  examined  in  details.  The  generation  of  electron  - 
positron  pairs  during  runaway  breakdown  will  be  studied  as  well.  It  would  be  demonstrated, 
that  the  last  process  is  quite  effective  what  allows  to  predict  the  possibility  to  observe  e^e” 
annihilation  y-line  during  powerfull  atmospheric  discharges. 


1.3  Kinetic  theory  of  runaway  breakdown  at  high  energies 

At  high  energies  to  describe  the  runaway  breakdown  it  is  necessary  in  addition  to 
ionization  effect  take  into  account  bremsstrahlung  processes  for  electrons,  additional 
ionization  by  high  energy  photons  and  Compton  effect  for  photons.  To  estimate  the 
influence  of  all  these  processes  let  us  first  discuss  kinetic  equations  for  electrons  and 
photons  in  a  simple  "x"  approximation 
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Here  rie  and  tiph  -  number  densities  of  electrons  and  photons,  Nm  and  Z  -  number  density  of 
molecules  and  full  charge  of  their  nucleus,  a  =  e  /he  is  a  fine  structure  constant,  ro=e  /me  - 
classical  radius  of  electron.  The  first  term  in  (1)  describes  the  absorbtion  of  photons  because 
of  Compton  effect,  4  is  Compton  absorption  length,  the  second  term  in  (1)  describes  the 
photon  generation  due  to  bremsstrahlung  (Ibr  is  the  length  of  generation  of  photon).  The  first 
term  in  (2)  describes  the  ionization  of  neutral  moleculas  by  fast  electrons  (  //  -  ionization 
length,  I(Z)=13.5  Z  eV),  and  the  second  one  by  photons. 

Let  us  first  compare  the  ionization  and  Compton  terms  (3).  Analyzing  the  equations  (1)  and 
(2)  it  is  possible  to  show,  that  in  conditions 


rln^(2r)  «  (rdaZ)  lr?[f^  me^/I(Z)] 


(4) 


the  avalanche  is  determined  mostly  by  ionization  process.  In  this  work  condition  (4)  is 
supposed  to  be  fulfilled  always. 

The  full  kinetic  equation  describing  the  runaway  breakdown  process  takes 
inconditions 

(4)  the  following  form  [6]. 
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Here  f(p,fAt)  is  the  electron  distribution  function,/?  -  electron  momentum,  E -  electric  field, 
p  =  (E  p}/(E  p)  -  cosine  of  pitch  angle.  The  integral  St  describes  collision  of  fast  electrons 
with  atomic  electrons  and  nucleus.  Fd  and  Fbr  is  friction  forces  of  electron 
determining  energy  losses  due  to  ionization  and  bremsstrahlung: 


\  at  j  me 

Here  Z,  is  nucleus  charge. 

The  second  term  in  (5)  describes  pitch  angle  scattering  of  runaway  electrons.  Term  Si 
describes  appearance  of  a  new  born  runaway  electrons,  y  -  Lorentz  factor  of  electron,  s  and 
W  -  kinetic  energy  of  fast  electron  s  =  mc^(y- 1 ),  W=  mc^(y'  -  1 ). 

At  Fig.  1.1  the  full  friction  force  is  demonstrated.  One  can  see,  that  the  ionization 
plays  the  main  role  at  low  energies  £•  <5  MeV.  At  £•  >  5  MeV  the  dominant  role  belongs  to 
bremsstrahlung  force.  The  basic  relation  of  the  theory  which  determines  the  runaway 
boundaries  eE  =  F  now  has  two  roots.  The  first  root  y  =  yi  determines  as  usual  critical 
electron  energy  Sc  [1].  The  runaway  region  lays  between  yj  and  y2  .  The  second  root  is 
determined  mainly  by  bremsstrahlung  force.  The  runaway  boundaries  energies  8c=mc^(yi-l) 
and  £c2=mc^(y2-l)  depend  on  relation  E/Ec 
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Fmin  is  the  minimal  drag  force  of  electrons  [1].  One  can  see  from  (7),  that  the  width  of 
runaway  region 


Tia 


me 
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The  highest  energy  £c2  which  can  be  reached  by  runaway  electron  is  growing  with 
E/Ec  approximately  linearly:  at  E/Ec  =5  it  reaches  the  value  430  MeV,  at  E/Ec=10  -  the 
value  850  MeV.  This  values  are  fully  determined  by  bremsstrahlung  force.  We  note,  that  in 
previous  papers  the  term  (6)  was  omitted.  It  means  that  kinetic  theory  was  considered  in  [1] 
for  low  energy  diapason  e<«ec2  only. 

Analyzing  collision  integral  we  see  that  pitch  angle  scattering  due  to  relativistic 
growth  of  electron  mass  with  energy  increase  falls  down  more  efficient  than  drag  force.  So, 
it  could  be  neglected  if  conditions 


^«1  (9) 

are  fulfilled.  At  such  conditions  pitch  angle  scattering  is  not  significant  and  due  to  this  the 
distribution  function  becomes  strongly  directed  along  electric  field  (see  [7]).  Thus  one  can 
consider  one  dimensional  -  integrated  over  angles  distribution  function. 
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Introducing  dimensionless  variables  yt. 
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Here  Se=E/Ec  is  electric  field  parameter,  function  <pd: 
determines  dimensionless  ionization  losses  and: 


mi 


-  bremsstrahlung  losses. 

The  ionization  term  takes  a  simple  form: 
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Thus  in  a  high  energy  range  due  to  relativistic  growth  of  the  electron  mass  its 
scattering  by  atomic  neucleus  is  weak  enough  and  the  distribution  function  of  energetic 
electrons  became  strongly  alongated  along  electric  field.  It  allows  to  use  kinetic  equation  in 
a  more  simple  form  (1 1),  (12). 


1.4  Time  development  of  distribution  function 

The  threshold  field  for  runaway  breakdown  Ec  is  usually  much  less  than  for 
conventional  breakdown  Et.  As  an  example  for  air  and  for  hydogenium  Ec/Eth  ~  0.1.  On  the 
other  hand  to  excite  runaway  breakdown  one  needs  fast  electrons  with  energies  s  >  Sc.  That 
determines  an  important  difference  between  runaway  and  conventional  breakdown.  How 
the  time  development  of  distribution  function  depends  on  the  existance  of  fast  electrons  or 
more  generally  on  initial  function  fo(yP.  The  problem  was  studied  numerically.  The 
numerical  soutions  were  obtained  in  runaway  region  (y  <  jci)  by  predictor  -  corrector  and 
Crank  -  Nickolson  methods  with  boundary  condition  f(yc^=0.  Both  methods  demonstrates 
identical  results.  To  determine  the  dependance  of  runaway  process  on  the  energy  of  seed 
electrons  the  initial  function  fo(y)  was  taken  usually  close  to  delta  function  fo(y)  ^S(y-f) . 

Results  of  numerical  calculations  of  equation  (1 1)  are  demonstrated  at  the  Fig.  1.2  - 
Fig.  1.5.  We  see  from  Fig.  1.2,  that  aXf<yi  (or  s*  <  s)  independing  on  the  value  of  electric 
field  {E  >  Ec),  the  initial  distribution  function  dissapears  in  time  and  initial  number  density 
of  electrons  Ne  falls  down  exponentially.  So,  without  seed  fast  electrons  runaway 
breakdown  does  not  exist. 

On  the  contrary,  if  f  >  yi  the  runaway  breakdown  is  effectively  developing.  Its 
main  features  seen  from  Fig.  1.3:  independently  on  the  value  f  the  distribution  function  at 
low  energies  develops  very  fast  to  its  quasistationary  state 

f(yt)=Fo(y)  (13) 

The  typical  time  of  this  process 

t,~X'  (14) 

At  the  high  energy  range  (y  >  y*)  the  establishing  of  quasistationary  solution  is  much  more 
slow  and  takes  time 

t2-(y2  -  y*)/(SE  - 1)  Se>1.2  (15) 
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Note  that  the  time  t2  is  diminishing  with  the  growth  of  y*  and  at  y*  >  y2  its  became  an  order 
of  ti. 

The  established  asympthotic  solution  takes  the  exponential  form  (13)  quite 
analogous  to  the  conventional  breakdown.  It  is  natural,  because  runaway  breakdown  is 
described  by  linear  equation.  Stationary  asympthotic  distribution  function  Fo(y)  is  shown  at 
Fig.  1.4.  for  different  values  of  parameter  Se.  One  can  see,  that  function  Fo  falls  down 
monotonically  with  y,  sharply  increasing  at  small  energies  s  -  mc^(y-l)  and  disappearing  at 
the  second  runaway  boundary  Fo(y2)  =  0.  The  behaviour  of  asympthotic  solution  Fo  could 
be  studied  analytically.  The  analytical  analysis  of  equation  (5)  shows  that  ionization  integral 
and  friction  term  are  most  essential  at  small  energies.  The  balance  of  these  terms  gives 
asymptotics  at  s  -^0 

In  intermediate  energy  range  yi  «  y  «  y2  the  distribution  function  falls  down 
exponentially 


/  exp' 


To  determine  cutoff  the  distribution  function  near  right  stationary  point 
y2  we  can  neglect  the  ionization  term.  Expanding  friction  force  near  this  point  we  can  find 
the  solution 


1  -—fjp  >0  p  =  mcy 

4 

These  results  are  in  agreement  with  numerical  simulation. 

1.5  The  avalanche  rate 

One  of  the  most  important  parameter  of  breakdown  is  the  avalanche  rate.  The 
increasing  of  number  density  of  electrons  N(t)  is  shown  for  different  beam  energies  y*  at 
Fig.  1.5.  We  see  that  the  asympthotic  solution  demonstrates  exponential  growth  with 
parameter  X: 

N(t)  =Noe^‘ 

We  investigated  the  dependence  of  avalanche  rate  X  on  dimensionless  electric  field  S 
.  The  results  are  shown  at  Fig.  1.6.  For  high  enough  values  of  ^this  behavior  is  well  fitted  by 
the  law: 

A  =  0.064  5^^^  (16) 

Note  that  as  seen  from  Fig.  1.5.  the  asympthotic  time  behaviour  of  electron  density 
does  not  depend  on  on  the  initial  distribution  function  (if  /*  >  y^  in  agreement  with  (13). 
Some  small  difference  in  }4(t)  is  seen  in  the  begining  of  the  process  at  time  scale  of  the 
order  of  ti  (14).  We  see,  that  the  full  time  of  establishment  of  distribution  function  (15)  is 
not  significant  —  electron  density  growth  is  going  with  the  same  avalanche  rate  from  the 
very  begining.  It  means  that  the  avalanche  process  practically  does  not  depend  on  the 
established  distribution  function  Foiy)  in  the  whole  breakdown  region.  The  dominant  role  is 
played  by  its  initial  part  in  the  vicinity  of  the  left  stationary  point  —  the  first  root  yi.  In 
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other  words,  only  the  region  near  critical  energy  Ec  determines  the  avalanche  rate.  This  fact 
is  in  agreement  with  qualitative  theory  [1]. 

Basing  on  this  analysis  we  can  estimate  analytically  the  avalanche  rate.  Really,  near 
left  stationary  point  it  is  possible  to  consider  nonrelativistic  approximation.  In  this 
approximation  the  kinetic  equation  (5)  is  reduced  to  the  form; 
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There  is  an  important  property  of  the  equation  (17)  -  it  allower  to  eliminate  the  electric  field 
Se.  Introducing  variables 


r  =  2tSp, 

one  can  find 
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Searching  exponential  solution  of  equation  {\%)fo=f  we  reduce  it  to  the  form 


(18) 
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Here  <I>  =  |f(w)wdw,  and  variable  y  =  =  u^S,,-. 

u 

Note,  that  parameter  a  determined  by  (5)  could  be  cosidered  in  equation  (19)  as  a  large 
parameter.  For  a  Doo  we  have  from  (19) 

(1_^)(P''_(I,'  =  0  <D  =  Coln(y-l)  +  C, 

This  solution  is  not  valid  for  ylJoo  and  yfJl .  To  match  this  solution  let  us  return  to  the  finite 
big  values  of  parameter  a  and  use  an  expanding  variables  introducing  z  as>’=  a^z. 

The  equation  (19)  takes  the  form: 


4(1  -  -  {cx4~z  + 1)0^  - =  0 

a  Aa  z 

For  aDoo  its  solution  is 
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Using  a  procedure  of  matching  solution  [8]  we  find  "combined"  solution  which  is  good  with 
exponential  accuracy  everywhere  except  the  vicinity  of  point  >>=7: 
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In  the  vicinity  of>’=/,  introducing  a:.- _y=  1+x  d'  we  find  the  solution 


<D  =  e  “ 

It  means  that  solution  (20)  is  cut  at  u*  =/+  4/a.  As  we  shall  see  below  this  cutting  defines 
the  avalanche  rate  with  logarithmic  accuracy.  Integrating  (18)  from  u*  to  oo,  one  can  find: 
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This  algebraic  equation  defines  parameter  X=2a8^'^/a.  For  a=ll  the  solution  of  (21)  gives 

X  «  0.055  5^'^ 

This  result  is  in  a  good  agreement  with  numerical  calculations  (16).  More  than  that  it  is  in 
agreement  with  previous  numerical  solution  of  kinetic  equation  in  the  low  energy  limit  [2-5] 
and  with  approximate  qualitative  theory  [1].  The  reason  of  this  agreement  is  in  the  fact  that 
main  role  in  ranaway  breakdown  process  plays  the  low  energy  runaway  boundary  £c. 


1.6  e*  e~  generation 

The  distribution  function  of  runaway  electrons  falls  down  rather  slow  (see(20)).  It 
means  that  fast  electrons  generates  intensive  flux  of  X  and  y  photons  due  to  bremsstrahlung 
process.  The  high  energy  photons  can  generate  electron  -  positron  pairs. 

The  process  of  X  and  y  generation  by  fast  electrons  in  neutral  gas  is  well  known. 
The  total  number  of  photons  could  be  estimated  using  emissivity  [9] 

®  ar^s 

at  £ 

Here  A/g  -  the  number  density  of  fast  electrons  which  generate  bremsstrahlung  emission  and 
-  density  of  neutrals. 

The  crossection  of  e*  generation  is  quite  analogous  to  Bethe  and  Heitler 
bremsstrahlung  crossection  [9] 

da  =  -  sin sin  edej(l)  X  D  (22) 

In  0)  q 

D  =  ^{4£l  -q^)s\nO^  +^{A£l  -q^)?Xn6_  — 6  )- 

Kl  K_  K^K_ 

-  (2g^  +2£'f^^)sin^^sin^_cos^ 
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Here 


a:±  =£’± -p±cos0^,  q^ -(p^+ p_-kY,  £^+£_=co 

Integrating  this  crossection  we  can  estimate  the  production 


dt 


=  <^±cNn,^r 


Taking  into  account  that  grows  exponentially  with  t,  we  find : 
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(23) 


This  relation  determines  the  number  of  generated  e‘^e“  pairs. 

To  estimate  the  emission  intensity  of  the  e^e~  annihilation  line  intensity  generated 
during  thunderstorm  it  is  natural  to  compare  with  the  results  of  observations  y  -  burst  from 
high  altitude  lightning  at  "Compton"  observatory.  According  to  Fishman  et  al  [10]  and 
Nemiroff  et  al.  [1 1]  the  flux  of  gamma  quants  in  energy  diapason  300  ^  1000  KeV  observed 
by  BATSE  at  the  heights  500  km  was  of  the  order  of  100  photons/cm^s.  Using  relation  (23) 
and  taking  into  account  that  N/Ne  ~  a~  10'^  one  can  estimate  the  flux  of  photons  from  e'^e” 
annihilation  at  the  distance  100  km  from  the  source  as  20  photons/cm^s.  Taking  into 
account  that  the  width  of  unihilation  line  is  narrow  one  can  suppose  that  the  y-quants 
generated  during  thunderstorm  could  be  detected  in  experiments. 
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Figure  1.1  Dependence  of  the  full  friction  force  F  on  y  (solid  curve).  The  dashed  curev 
presents  the  friction  force  without  accounting  bremsstrahlung. 
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2.  Investigation  of  runaway  breakdown  phenomena  in  laboratory  installation  in  conditions  of  cyclotron 
resonance 


Fig.2.1.  Scheme  of  the  experimental  setup:  1  -  magnetic  system;  2  -  vacuum  chamber;  3  - 
microwave  generator;  4  -  attenuator  or  isolator;  5  -  obturator;  6  -  waveguide 
components;  7  -  horn  coupler;  8  -  vacuum  pump;  9  -  water  attenuator  section;  10  - 
horn  antenna;  1 1  -  detector;  12  -  loop  antenna;  13  -  microwave  receiver  with 
spectrum  analyzer;  14  -  X-ray  receiver;  15  -  optic  fiber;  16  -  photomultiplier;  17  - 
data  display  system  (computers). 

The  insertion  demonstrates  the  magnetic  field  distribution  in  stationary  state 
normalized  to  magnetic  field  in  the  center  of  vacuum  chamber. 
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2. 1  Setup  description 

The  setup  represented  in  Fig.2.1  includes  a  vacuum  chamber  2  of  quasispherical 
form,  which  has  a  the  distance  between  a  pair  of  opposite  flanges  450  mm.  The  chamber  is 
placed  between  two  similar  coaxial  copper  coils  1  of  a  static  magnetic  field.  The  magnetic 
field  force  line  distribution  inside  the  chamber  has  configuration  of  a  bottle  with  two 
symmetric  necks  which  forms  a  magnetic  trap  for  collisionless  electrons.  The  trapped 
electrons  confinement  time  in  the  chamber  is  likely  to  be  enough  to  confine  a  bulk  of  initial 
electrons  of  large  energies  during  the  time  of  run-away  electrons  breakdown  in  air.  The 
magnitude  of  the  static  magnetic  field  at  the  center  of  chamber  varies  relative  to  value  875 
Gs  which  corresponds  to  the  ECR  for  the  electric  field  frequency  2,45  GHz. 

A  magnetron  3  generates  electromagnetic  oscillations  at  the  fixed  frequency  2,45 
GHz,  with  varying  power  up  to  5  kW  which  creates  the  ECR  accelerating  field  for  electrons 
caught  inside  the  mirror  trap.  The  microwave  power  is  transmitted  into  the  chamber  through 
a  waveguide  system  6  which  consists  of  the  adapter  with  pyrex  window  insulating  vacuum 
of  the  chamber  from  athmosphere,  the  horn  antenna  7,  the  circular  waveguide  with  travelling 
TEii  mode,  the  smooth  waveguide  transition  from  a  rectangular  cross  section  to  circular  one, 
the  rectangular  waveguide  elbow,  the  water  attenuator  or  ferrite  circulator  4.  For  purpose  of 
pulse  periodic  regime  realization  it  was  made  mechanical  waveguide  cutter  off.  The 
maximum  microwave  electric  field  strength  at  the  center  of  vacuum  chamber  under 
condition  of  travelling  wave  at  maximum  power  is  Eo«150  V/cm.  In  the  vicinity  of  the  ECR 
the  microwave  electric  field  is  equivalent  to  a  static  electric  accelerating  field  for  electrons. 
Following  experimental  techniques  for  study  of  x-ray  bremmsstrahlung  radiation  and 
optical  radiation  produced  by  the  run  away  electrons  and  microwave  synchrotron  radiation 
at  the  harmonics  of  the  accelerating  field  were  used: 

•  Scintillation  x-ray  transducer  14  (BYCRON,  model  2.5M2.5/3LP-X). 

The  transducer  consists  of  NaJ(Tl)  scintillator  covered  with  aluminium  cup  of  1  mm 
thickness  and  combined  with  a  photomultiplier.  The  optical  axis  of  the  transducer  placed  in 
the  chamber  was  perpendicularly  directed  to  the  trap  axis.  The  operating  quantum  energy 
range  of  the  x-  ray  transducer  100  keV  -h  1  MeV.  Its  signal  was  registered  at  a  monitor  of 
the  PC  «Pentium»  by  means  of  an  oscillograph  card. 

•  Photomultiplier  16  for  observing  of  visible  air  glow  discharge  under  action  of 
both  ECR-discharge  and  run  away  electrons. 

•  Personal  x-ray  doze  meter  (not  shown)  for  measurement  of  x-ray  spectrum 
outside  the  chamber.  The  operating  quantum  energy  range  of  the  counter  0,05  h- 
1,5  MeV.  Its  error  of  doze  power  measurements  is  +  30%;  exposition  time  is  36  s; 
the  dependence  of  doze  measurements  over  the  operating  energy  range  is  no  more 
45%. 

•  Microwave  spectrum  analyser  13  over  the  frequency  range  of  2-nd  and  3-d 
harmonics  for  the  base  frequency  of  the  microwave  generator.  Signals  for  the 
analyser  were  received  from  loop  antenna  12  at  end  of  a  coaxial  feeder 
introduced  along  the  trap  axis. 

•  Detector  1 1  for  monitoring  incident  MW  power. 

2.2  Experimental  Results 

In  the  experiments  with  the  described  setup  the  intense  X-ray  emission  was  found  to  be 
generated  under  ECR  discharge  at  the  pressure  range  <10"^  Torr  due  to  bremmsstrahlung 
radiation  of  electrons  accelerated  in  relatively  weak  electric  field  E=5-^20  V/cm  up  to 
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subrelativistic  energies  100-^300  keV.  Simultaneously  the  MW  synchrotron  radiation  at  the 
second  harmonic  was  observed.  Its  spectrum  points  out  the  presence  bulk  of  electrons  with 
energies  <10  keV  in  the  system.  Specific  feature  was  that  the  beginning  of  discharge  visible 
emission  owing  to  low  energy  part  formation  of  electron  spectrum  (<100  eV)  had  a  delay 
reaching  80  ms  in  respect  to  the  appearance  of  X-ray  and  synchrotron  signals.  This  allows 
to  separate  the  whole  process  of  the  ECR  breakdown  into  following  two  stages:  first  of 
which  to  be  an  initial  stage  of  breakdown  with  participating  of  high  energy  run  away 
electrons,  and  second  one  to  be  ordinary  ECR  discharge  with  participating  of  low  energy 
electrons.  Such  separation  is  reasonable  due  to  difference  of  physical  processes  resulting  in 
the  characteristic  of  electron  distribution  functions.lt  allows  to  find  analogy  of  ECR 
breakdown  with  the  model  of  runaway  electron  breakdown  by  Gurevich,  Milikh  and 
Roussel-Dupre  applied  for  explanation  of  lightning  mechanism  in  Earth  atmosphere. 
Recently  interesting  results  were  also  observed  in  pulse  mode  of  microwave  power.  Such 
mode  can  be  achieved  by  regular  cutting  off  the  incident  MW  power  with  mechanical 
waveguide  modulator.  The  preliminary  results  are  shown  in  Fig.2.2. 
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One  can  see  (Fig.2.2)  that  at  air  pressure  S  Torr  we  always  have  a  time  delay 
visible  light  radiation  of  discharge  in  respect  to  the  X-ray  signal  and  the  cut  off  transmitted 
MW  radiation  connected  with  that  process  of  accumulation  of  low  energy  plasma  with 
overcritical  electron  density. 

The  behavior  of  signals  in  Fig.2.3  corresponding  to  air  pressure  2T0‘^  demonstrates 
the  visible  discharge  is  ahead  in  respect  to  the  X-ray  signal.  However  one  can  see  the 
flickering  of  discharge  when  the  incident  microwave  energy  is  redistributed  between 
runaway  electrons  and  low  energy  plasma. 

The  main  objectives  of  our  future  work  are  detailed  investigations  of  runaway 
electron  breakdown  taking  into  account  its  temporal  and  spatial  dynamics,  definition  of 
energy  characteristics  of  the  process  and  energy  balance,  measurements  of  energy  y-quanta 
spectra  depending  on  experimental  conditions. 

To  study  the  dynamics  of  the  discharge  plasma  density  due  to  the  RB  a  microwave 
interferometer  (^=8  mm)  of  opened  resonator  type  will  be  constructed.  The  x-ray 
spectrometer  in  the  energy  range  30  keV  -  1  MeV  will  be  created  with  use  of  the  X-ray 
scintillation  detector  (in  our  possession)  and  the  pACE^^-system  (EG&G  ORTEC)  which 
consist  of  amplifier,  multichannel  analyzer  plug-in  card  with  emulation  software(will  be 
ordered  to  the  ISTC).  The  averaged  x-ray  spectra  depending  on  the  regimes  of  the  RB  and 
optical  spectra  of  the  RB  plasma  will  be  studied.  Both  probe  measurements  and  electron 
energy  distribution  function  measurements  with  an  electrostatic  analizer  for  electrons 
leaving  the  magnetic  trap  will  be  carried  out. 

We  intend  to  continue  investigation  with  pulse  mode  of  MW  generation  because  it’s 
reasonable  to  have  the  time  limited  action  on  the  RB  for  study  of  its  dynamics.  For  this 
purpose  the  waveguide  modulator  5  (Fig.2.1)  cutting  off  continuous  microwave  radiation 
was  constructed.  It  allows  to  vary  duration  and  period  of  pulses  in  wide  range.  The  demands 
for  the  modulator  was  caused  by  the  reason  that  MW  generator  is  able  to  work  in 
continuous  mode  generation,  but  not  able  to  do  it  in  pulse  mode  if  to  use  pulse  high  voltage 
supply.  Microwave  pulses  after  the  modulator  have  trapezoidal  form  with  rise  and  fall  time 
of  order  1/3  of  plateau.  Owing  to  the  pulse  mode  we  hope  to  study  lifetime  of  high  energy 
electrons  caught  by  the  magnetic  trap  during  a  previous  microwave  pulse  and  their  role  as 
seed  in  the  RB  phenomenon  during  followed  pulse. 
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3.  Investigation  of  an  influence  of  combined  effects  of  runaway  breakdown  and 
cosmic  rays  on  lightning  processes  in  thunderstorm  atmosphere 

3. 1  Abstract 

The  intensive  X-ray  emission  during  thunderstorms  was  studied  at  the  height  3340  m  above 
sea  level.  For  the  first  time  the  short-time  (1 — 5  min)  bursts  of  X-ray  emission  were 
observed.  The  bursts  are  highly  correlated  at  the  wide  space  region  (about  0.5  km).  The 
main  component  of  X-ray  emission  in  bursts  has  the  energy  50 — 80  keV.  The  observed 
bursts  could  be  attributed  to  the  bremsstrahlung  determinated  by  the  runaway  breakdown 
(RB)  effect  in  the  thundercloud  electric  field.  Ground  based  observations  of  RB  open  a  wide 
opportunities  for  the  studies  of  fundamental  processes  in  thunderstorms. 

3.2  Introduction 

The  new  mechanism  of  lightning  initiation  in  thunderstorm  clouds  based  on  a 
combined  action  of  runaway  breakdown  (RB)  and  extensive  atmospheric  showers  (EAS) 
was  proposed  in  [1].  According  to  this  theory  in  shower  core  the  ionization  of  the  air  due  to 
RB  effect  is  growing  strong  enough  to  produce  a  local  highly  conductive  plasma.  The 
thunderstorm  electric  field  is  amplified  strongly  inside  such  region.  That  can  serve  to  the 
triggering  of  a  lightning  discharge. 

The  Lebedev  Physics  Institute  (LPI)  group  proposed  an  experimental  program  to 
study  this  mechanism.  The  program  includes  simultaneous  observations  of  cosmic  ray 
showers,  runaway  breakdown  emission  and  lightning  discharges  at  the  mountain  cosmic  ray 
station  (MCRS).  The  LPI  MCRS  is  situated  on  a  pass  in  Tien-Shan  mountains  at  the  height 
3340  m  above  sea  level  (effective  atmosphere  depth  690  g/cm^).  The  thunderstorm  clouds 
go  through  the  pass  at  the  height  about  0.1 — 1  km  during  of  three  months  yearly.  The 
installation  can  detect  extensive  air  showers  measuring  simultaneously  EAS  electrons  with 
energies  0.1-^10  MeV,  muons  with  energies  above  5  GeV  and  neutrons  with  energies  above 
100  MeV.  The  RB  effect  lead  to  a  strong  amplification  of  the  X-ray  emission  inside  and 
below  the  thunderclouds  [2]-[6].  To  observe  RB  emission  a  system  of  Geiger  counters  was 
installed  in  a  wide  region  (approximately  500  m  around  MCRS).  To  fix  the  position  of 
lightning  discharge  by  its  strong  electromagnetic  and  acoustic  radiation  a  group  of  special 
antennas  together  with  directed  microphones  was  established. 

In  the  present  work  we  inform  about  the  first  results  of  X-ray  emission 
measurements.  The  strong  amplification  of  emission  during  the  whole  thunderstorm  period 
is  found.  For  the  first  time  the  bursts  —  short  time  strong  peaks  of  X-ray  radiation  were 
observed.  The  typical  time  length  of  such  a  burst  is  l-;-5  min  and  radiation  increases  are 
highly  correlated  in  different  points  at  the  scales  minimum  of  the  order  of  several  hundred 
meters. 

3.3  Experimental  Installation 

The  schematic  layout  of  MCRS ’s  detectors  disposition  is  shown  at  Fig.3.1.  The 
installation  used  in  this  experiment  consists  of  the  following  main  parts: 

•  The  shower  system  on  the  base  of  scintillation  detectors  located  nearly  concentric 
circumferences  on  the  flat  part  of  the  pass  (the  “Hadron”  experiment  installation  [7]). 
The  scintillators  can  detect  the  high  energy  (above  10  MeV)  electrons  and  gamma  quanta 
travelling  through  the  installation  within  an  extensive  air  shower  fronts. 

•  The  NM64  type  neutron  supermonitor  [8]  consisting  of  three  independant  units.  Each 
unit  has  3  m^  sensitive  area  and  includes  six  registration  channels  on  the  base  of  neutron 
counters  SNM15.  The  supermonitor  was  designed  for  the  registration  of  one-minute 
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intensity  variations  of  cosmic  ray  hadrons  with  energy  above  50 — 100  MeV.  The 
electronics  of  neutron  monitor’s  registration  channels  is  designed  to  provide  methodic 
error  at  the  level  of  0.01%  in  intensity  measurements.  The  input  threshold  value  for 
neutron  counter's  pulses  is  set  to  be  equal  to  10  mV,  corresponding  background  pulse 
intensity  being  2700  pulses  per  min  (p.p.m.)  for  each  unit. 

•  The  boxies  of  high  threshould  gamma  ray  and  electron  detector  are  placed  above 
supermonitor  units  as  shown  at  Fig.S.lb.  This  detector  is  made  on  the  base  of  Geiger- 
Mueller  counter  tubes  of  SI5G  type  which  are  capable  to  register  hard  X-  and  soft  y-ray 
emissions  with  nearly  1%  efficiensy.  Four  detector  boxies  each  of  wich  contains  20 
counters  with  parallel  outputs  and  corresponds  to  the  one  registration  channel  were 
mounted  above  each  of  three  supermonitir  units,  their  total  sensetive  area  being  equal  to 
8  m^.  Because  the  counters  of  this  detector  are  operating  in  proportional  regime  it  is 
necessary  to  decrease  their  dead  times  up  to  5-10  mcs  (which  is  20  times  shorter  as  for 
usual  Geiger  counters).  To  achive  the  charged  particles  registration  efficiency  as  high  as 
for  Geiger  regime  (above  95%)  the  registration  channel's  threshould  is  set  to  20  mV. 
The  roof  of  the  building  where  neutron  supermonitor  and  high  threshould  y-ray  detector 
is  situated  is  made  of  0.8  mm  thick  iron,  the  wood  ceiling  has  thikness  about  15  cm  and 
the  detector’s  boxes  are  made  of  2  mm  thick  aluminium.  It  provides  the  upper  threshold 
about  25 — 30  MeV  for  relativistic  charged  particles  registration.  The  main  absorber 
above  the  detector  is  carbon  with  the  effective  thikness  13.6  g/cm^.  Thus,  the  threshold 
value  for  the  photons  is  about  30  keV.  The  general  count  intensity  for  all  channels  is 
determinated  mainly  by  the  electron  background  flux  and  is  equal  to  3.210^  p.p.m.  The 
corresponding  statistical  error  in  1 -minute  data  is  less  than  0.2%. 

•  The  muon  detector  made  of  SI5G  type  counters  with  common  sensitive  area  2.7  m^  is 
placed  under  one  of  the  supermonitor  units  (see  Fig.3.1b).  Output  pulses  from  each  of 
four  lower  detector  boxies  are  switched  to  a  coincidence  scheme  with  the  pulses  coming 
from  corresponding  boxies,  located  above  supermonitor  unit.  The  thickness  of  matter 
above  lower  boxies  is  15  g/cm^  of  carbon  and  150  g/cm^  of  lead  absorber.  The  latter 
exeeds  35  electromagnetic  cascade  units,  excluding  completly  the  registration  of 
atmospheric  photons  and  electrons  with  energies  below  10  GeV.  The  detector  registrates 
mainly  the  muons  with  energies  above  220  MeV  travelling  through  the  unit  in  nearly 
vertical  direction.The  general  count  intensity  is  about  7-10'^  p.p.s.  Hence,  the  fluctuations 
in  one-minute  data  are  about  0.5%.  The  registration  channel’s  characteristics  of  muon 
detector  are  just  the  same  as  those  for  the  high  threshould  gamma  detector. 

•  The  newly  build  low  threshould  wide  spread  X-ray  detector  (WS)  created  in  Summer 
1999  is  designed  to  investigate  the  thunderstorm  clouds  low  energy  electromagnetic 
emission.  It  consists  of  four  registration  points  situated  perpendicular  to  the  preferential 
direction  of  thunderstorm  clouds  motion  over  the  cosmic  ray  station's  pass  (see  Fig.3.1a). 
The  distance  between  the  extreme  right  and  extreme  left  points  is  about  500  m.  Each  of 
the  points  has  two  layers  of  sixty  SI5G  counters.  The  outputs  of  the  counters  in  each 
layer  are  combined  in  one  registration  channel.  Hence,  there  are  two  registration 
channels  in  each  point  with  sensitive  square  2  m^  each.  The  SI5G  counters  are  put  in 
Geiger-Mueller  regime,  and  the  threshold  of  registration  channels  is  equal  to  2 — 3  V. 
Electronic  circuits  in  each  detector  point  are  based  on  vacuum  tubes  in  order  to  protect 
them  from  the  nearby  lightning  discharges.  The  powerfull  vacuum  tubes  ensure  the 
output  pulse  front  length  not  worse  then  2 — 3  mcs  at  output  resistance  of  240  Ohm  and 
capacity  about  20  mcF.  Due  to  this  distortion  of  the  pulse  form  during  its  passage 
through  300 — 400  m  long  cable  is  negligible.  The  detectors  of  each  registration  points 
are  placed  within  a  little  light  house.  The  counters  are  separated  from  the  atmosphere  by 
the  0.8  mm  (or  0.74  g/cm^)  thick  iron  roof 
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All  the  detectors  described  above  have  the  bouth  independed  power  supply  and  data 
recording  systems.  Distant  detectors  are  connected  to  a  central  registration  system  (located 
in  the  building  with  neutron  supermonitor)  by  screened  underground  cables.  The  data  are 
transfered  to  computer  by  the  means  of  a  CAMAC  system.  To  avoid  the  danger  of 
semiconductor  electronics  damage  in  the  moment  of  lightning  discharge,  the  cables  are 
conjugated  with  CAMAC  scaler  blocks  throughout  vacuum  tube  based  circuits. 

The  pulses  coming  from  all  of  eight  channels  of  low  threshould  WS  X-ray  detector 
were  counted  by  CAMAC  scalers  during  of  one  second  time  intervals  and  accumulated 
numbers  of  them  were  recorded  by  registration  computer.  Amounts  of  pulses  obtained  from 
other  detectors  were  counted  during  of  one  minute  intervals  and  registrated  together  with 
the  clock  indication,  temperature  and  barometric  pressure.  The  registration  was  perfomed 
continuously  during  of  4  monthes  in  July — October  1999. 

We  have  tested  the  electric  noise  in  all  channels  by  switching  off  the  high  voltage 
power  supply  in  registration  points.  In  such  a  case  the  pulses  from  the  channels  were  not 
registrated  even  in  thunderstom  conditions. 

3.4  Experimental  results 

Intensity  variations.  During  of  July  and  up  to  the  last  days  of  August,  1 999  there 
were  no  thunderstorms  in  the  vicinity  of  mountain  station.  The  typical  measurements  of 
pulse  intensity,  obtained  from  various  detectors  at  this  quiet  period  are  presented  at  the  3.2. 
The  shown  data  correspond  to  the  total  intensity  registrated  by  all  channels  of  each  detector 
(i.e.  the  amounts  of  pulses  obtained  from  8  channels  of  low  threshould  X-ray  detector,  12 
channels  of  high  threshould  y-ray  detector,  3  neutron  monitor  units,  4  muon  coincidence 
channels  show  the  intensity  curves  at  this  figure). 

One  can  see  that  both  low  and  high  threshould  y-ray  detectors  exhibit  no  significant 
changes  in  count  velocity  through  the  whole  period  shown  (6  houres).  The  signals  from 
neutron  monitor  and  muon  detector  demonstrate  the  same  stability  (their  slow  trend  is 
caused  by  variation  of  atmospheric  pressure). 

In  the  end  of  August  1999  we  have  observed  five  cases  when  at  low  heights  over  the 
MCRS  were  passing  the  thunderstorm  cloud  fronts  accompanied  by  the  strong  lightning 
discharges.  The  same  data  as  given  at  Fig.3.2  for  a  quiet  day  are  presented  at  Fig.3.3a-e  for 
typical  thunderstorms  events.  Time  scale  at  the  bottom  of  the  figure  is  centered  relatively  to 
the  moment  when  thunderstorm  beginning  was  fixed.  This  moment  is  marked  by  a  vertical 
dotted  line.  Below  we  shall  discuss  the  typical  thunderstorm  event  on  August,  31. 
Corresponding  data  are  presented  at  Fig.  3.3.e. 

It  can  be  seen  from  the  typical  curves  shown  at  Fig.3.3e,  that  each  time  after  the 
moment  of  thunderstorm  beginning  a  dramatic  growths  of  emission  intensity  measured  both 
by  low  and  high  threshould  y-ray  detectors  occures.  We  emphasize  the  significant  difference 
between  these  two  types  of  detectors;  low  threshold  detectors  measured  practically  free  flux 
of  X-rays,  and  high  threshold  detector  that  was  under  the  roo/having  additional  absorption 
equivalent  to  13.6  g/cm^  of  carbon.  One  can  see  from  the  Fig.3.3.e,  that  the  signals  of  the 
free  (low  threshould)  and  under  the  roof  (high  threshould)  detectors  are  well  correlated  in 
time.  These  correlations  are  quite  regular  having  the  same  character  for  every  thunderstorm. 
The  characteristic  scale  of  long  time  emission  increses  is  2 — 3  houres  after  the  moment 
when  thunderstorm  clouds  begin  to  pass  the  MCRS. 

At  the  same  time  the  neutron  monitor  and  muon  deteetor  (which  are  sensitive  to 
energetic  particles  of  cosmic  radiation  only)  do  not  demonstrate  any  thunderstorm  caused 
effect. 

The  results  obtained  for  the  all  of  five  observed  thunderstorms  are  quite  analogous. 
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Space  correlations.  The  y-xay  intensity  increases  during  thunderstroms  discussed 
above  have  been  studied  in  data  streams  obtained  separately  from  each  of  low  threshould 
detector  points  shown  at  Fig.3.1a.  The  correlation  coefficients  between  all  available  pairs  of 
detector  points  were  determined.  The  results  —  time  dependences  of  correlation  coefficients 
Cjj  calculated  for  100  sec.  Length  time  periods  for  three  detector  points  pairs  —  are 
presented  at  Fig. 3.4  for  one  of  the  observed  thunderstorms  (the  corresponding  data  obtained 
for  all  other  detector  pairs  and  of  other  thunderstorms  are  quite  analogous).  One  can  see 
from  the  Fig.3.4  that  under  the  fine  weather  conditions  the  data  at  different  observational 
points  are  practically  noncorrelated.  During  the  thunderstorm  period  the  picture  is 
absolutely  different:  correlation  coefficients  for  all  low  threshold  detector  points  pairs 
becomes  very  high  —  about  0.8  and  even  higher.  This  high  correlation  demonstrate  that 
the  observed  of  X-ray  intensity  strong  growth  is  deeply  connected  with  thunderstorm 
processes  developing  simultaneously  at  the  wide  space  region  of  the  order  of  0.5  km,  maybe 
even  larger. 

Short  time  emission  bursts.  The  accompanying  thunderstorm  emission  intensity 
variations  are  well  correlated  in  time  also.  It  can  be  seen  from  Fig.3.3.e  for  X-ray  signals  of 
the  free  (curve  1)  and  under  the  roof  (curve  2)  detectors.  Analogous  time  correlations  exist 
for  the  emission  at  different  low  threshold  detector  points.  These  correlations  are  quite 
regular,  they  have  the  same  character  for  every  thunderstorm.  The  characteristic  timescale 
of  this  longtime  emission  variation  is  30-h60  min.  One  can  mention  the  following  peculiarity 
of  emission  intensity  —  it  has  often  a  double  or  triple  pike  structure  (see  Fig.3.3.e). 

An  interesting  new  feature  is  short  time  bursts  of  X-ray  emission.  Characteristic 
length  of  the  bursts  is  1-^5  min.  The  typical  short  time  bursts  occured  during  of 
thunderstorm  event  of  31  August  1999  can  be  seen  at  Fig.3.3.e  in  time  moments  0.2  (peak 
1),  0.4  (peak  2),  1.4  (peak  3)  houres  after  the  beginning  of  thunderstorm.  It  is  obvious  good 
correlation  between  the  intensity  bursts  of  low  and  high  threshould  emission  detectors.  Any 
increases  at  corresponding  time  moments  in  the  intensity  of  energetic  cosmic  ray 
components  (i.e.  neutrons  and  muons)  are  absent. 

The  one  of  the  bursts  data  for  low  threshould  detector  intensities  (peak  3)  are  shown 
in  large  scale  at  Fig.3.5  both  for  the  sum  of  all  channels  and  for  each  detector  point 
separately.  One  can  conclude  from  this  figure,  that  emission  intensity  bursts  are  correlated 
strongly  between  various  detectors  separated  with  space  distancies  of  the  order  of  some 
hundreds  of  meters.  According  to  statistical  errors  marked  at  Fig.3.5,  the  value  of  intensity 
increase  during  of  a  typical  emission  burst  exeeds  the  level  of  statistical  flluctuations  by 
15 — 30  times. 

Absorption  of  emission.  The  two  types  of  measurements  of  emission  intensity  were 
performed  simultaneously:  low  threshold  and  high  threshold,  which  demonstrate  quite  a 
good  time  correlation  (Fig.3.3.e).  The  following  main  difference  between  these 
measurements  exists.  The  low  threshold  detectors  were  separated  from  the  atmosphere  by 
the  thin  (0.8  mm)  iron  protector  only.  The  high  threshold  detectors  had  the  same  protector, 
but  they  were  installed  inside  the  neutron  supermonitor  building,  which  has  a  roof  with 
effective  thickness  13.6  g/cm^  of  carbon.  The  roof  gave  a  significant  additional  absorption. 
By  the  comparison  of  the  data  from  free  (low  threshold)  and  under  the  roof  (high  threshold) 
detectors  the  characteristic  absorption  of  emission  could  be  found. 

To  perform  comparison  we  determined  first  five-minute  averaged  background  of 
both  sets  of  data.  As  fast  electrons  determine  the  major  part  of  the  counts  we  have  to 
eliminate  the  background  and  compare  the  excess,  which  is  due  to  X-ray  effect.  The  excess 
relation  r  for  both  sets  of  data  is  presented  at  Fig.3.6.  One  can  see,  that  the  typical  value  of 
excess  relation  of  the  free  (low  threshold)  to  the  under  roof  (high  threshold)  data  is  2h-4. 
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Such  an  absorption  in  13.6  g/cm^  carbon  is  characteristic  for  energy  of  y-quanta  200-5-500 
keV[18]. 

It  is  important,  that  during  short  time  bursts  the  value  of  r  is  always  significantly 
higher,  growing  up  to  8-5-12  (Fig.3.6).  That  shows  that  the  energy  spectrum  of  X-rays  during 
short  bursts  is  essentially  different  from  the  main  part  of  thunderstorm:  its  major  component 
has  low  energy  X-rays  (<  50-5-80keV)  which  are  effectively  absorbed  by  the  roof  Note  that 
peaks  are  more  pronounced  in  absorption  curves  (Fig.3.6),  than  in  intensity  (Fig.3.3.e).  It 
opens  a  possibility  to  use  absorption  as  an  effective  method  of  X-ray  burst  detection. 

Flux  of  X-ray  quanta.  To  determine  the  flux  of  X-ray  quanta  we  have  first  to  average 
the  excess  counts  of  low  threshold  data  and  then  takes  into  account  absorption  by  the  iron 
protector  and  sensitivity  of  Geiger-Muller  counters  SI5G  for  X-ray  quanta.  That  gave  an 
approximate  value  of  flux  in  bursts  10-5-20  photons/(cm^  sec  sr). 

3.5  Discussion 

The  aim  of  this  work  is  to  study  the  influence  of  thunderstorms  on  the  X-ray  emission  from 
atmosphere.  Effective  growth  of  the  emission  in  all  observed  thunderstorms  is  definitely 
established. 

In  fine  weather  conditions  weak  disturbances  of  the  emission  (~  1%)  determined  partly  by 
the  changes  in  atmospheric  pressure  are  regularly  observed  (Fig.3.1)  That  agrees  well  with 
the  previous  statements  [9]  -  [11].  During  thunderstorm  the  picture  changes  fundamentally. 
The  X-ray  emission  is  strongly  amplified  for  hours.  The  emission  amplification  in  the 
thunderstorm  period  was  studied  previously  in  a  number  of  papers  (see  [10]-[12]). 

The  main  peculiarities  of  our  observations  are: 

1 .  X-ray  emission  was  studied  at  the  ground,  but  high  in  the  mountains  (height  3340m). 
Analogous  investigations  were  done  previously  only  by  Shaw  [10]  and  Suszcynsky  et 
al.  [1 1]  (mainly  at  lower  heights). 

2.  For  the  first  time  the  emission  was  observed  simultaneously  in  several  points  with  the 
large  distance  between  them  (Fig.3.1). 

3.  For  the  first  time  the  free  and  partly  screened  emission  was  studied  simultaneously. 

Below  we'll  compare  our  observations  with  the  previous  studies.  The  new  results  will  be 
discussed. 

The  main  mechanisms  of  strong  amplification  of  X-ray  emission  intensity  during 
thunderstorm  are: 

1 .  Gamma  ray  emission  from  the  beta  decay  of  radon  daughter-ion  products. 

2.  Bremsstrahlung  of  accelerated  electrons  due  to  the  action  of  a  strong  thunderstorm 
electric  field. 

Radon  effect  is  determined  by  the  precipitation  of  radon  and  its  products  to  the  ground 
during  rainfall.  Gamma  ray  emission  from  radon  daughter-ion  [3  decay  has  energies  of 
quanta  0.2-5-2  MeV.  The  decay  processes  and  rainfall  determine  its  main  time 
characteristics:  20  -5-  30  min  growing  up  and  about  50  min  after  washout  [1 1,13]. 

Electric  field  effect  is  determined  by  the  acceleration  of  electrons.  Wilson  [14]  was  the  first 
who  suggested  that  the  secondary  electrons  generated  by  cosmic  rays  could  be  accelerated 
by  the  thunderstorm  electric  field.  A  new  fundamental  step  was  done  in  the  runaway 
breakdown  (RB)  theory  developed  by  Gurevich  et  al.  [2].  It  was  shown  that  accelerated 
electrons  could  organize  an  exponentially  growing  avalanche  in  thunderstorm  electric  field. 
RB  effect  is  supposed  to  determine  maximal  electric  field  reached  in  thunderclouds 
[2,15,16].  Directly  RB  process  we  discussed  in  the  Introduction  in  connection  with  the 
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lightning  triggering  theory  [1].  RB  effect  lead  to  the  strong  amplification  of  X-ray  emission 
due  to  the  bremsstrahlung. 

Amplification  of  X-ray  emission  attributed  to  RB  effect  was  observed  in  [3,4,17]  during  the 
measurements  inside  thunderclouds  at  the  balloons.  The  main  features  of  these  observations 
could  be  summarized  as  follows: 

a)  The  intensity  of  X-ray  emission  at  the  height  ~  4  km  is  sharply  (in  seconds  time  scale) 
growing  up  to  10^  in  the  comparison  with  the  standard  level  determined  by  cosmic  rays. 
The  maximum  of  observed  emission  bursts  lays  in  the  diapason  60-^90  keV,  the 
maximal  count  rate  is  40^50  photon/(cm^  sec  sr)  [4,17]. 

b)  The  horizontal  scale  of  the  region  of  emission  amplification  is  of  the  order  1-km  or 
more  in  agreement  with  the  layered  structure  of  thunderstorm  clouds. 

c)  The  time  scale  of  X-ray  bursts  is  of  the  order  1  min. 

d)  The  intensity  of  bursts  diminishes  sharply  (in  second  time  scale)  with  the  lightning.  But 
in  a  10-sec  time  scale  it  can  grow  up  again. 

e)  Lightning  was  seen  during  all  bursts.  So,  it  is  natural  to  suppose,  that  a  correlation  of 
lightning  events  with  X-ray  intensity  bursts  do  exist. 

Let  us  discuss  now  results  of  our  observations,  taking  into  account  the  main  features  of 
both  considered  mechanisms  of  X-ray  emission  amplification. 

1.  The  general  structure  of  a  long  time  emission  (Fig.3.3.e)  agrees  well  with  the  radon 
source.  A  peculiarity,  which  could  be  mentioned,  is  the  usual  double  pike  structure  of 
our  intensity  curves  (Fig.3.3.e).  The  energy  of  gamma  quanta  (0.2h-0.5  MeV), 
determined  by  the  absorption  in  the  13.6  g/cm^  carbon  for  the  main  part  of  observed 
thunderstorm  emission  (Fig.3.6),  agrees  with  radon  source  also. 

2.  Short  time  bursts  of  X-ray  emission  (Fig.3.3.e-3.6)  are  the  main  new  result  obtained  in 
our  measurements.  The  bursts  are  well  correlated  at  the  scales  0.5-km  or  more.  Such 
correlated  short  time  bursts  in  the  ground-based  measurements  are  observed  for  the  first 
time.  There  exists  a  quite  good  correlation  between  the  observations  of  free  and  partially 
absorbed  bursts.  The  analysis  of  absorption  effect  shows,  that  the  energy  spectrum  of  X- 
rays  in  the  bursts  is  essentially  different  from  the  main  part  of  thunderstorm:  its  major 
part  consists  of  the  low  energy  X-rays  (E  <  50-^80  keV).  Such  type  of  energy  spectrum 
is  characteristic  for  RB  effect  only. 

In  Table  1  the  main  features  of  our  observations  of  bursts  are  compared  with  the  X-ray 
emission  measured  at  balloons.  We  see  that  the  features  of  the  balloon  measurements 
attributed  to  RB  effect  are  rather  close  to  the  short  burst  observed  on  the  ground  in  our 
measurements  (Table  1).  One  can  conclude  that  in  short  bursts  the  runaway  breakdown 
bremsstrahlung  was  observed. 

We  note  that  this  statement  does  not  contradict  to  the  other  ground  based  studies,  including 
the  mostly  comprehensive  and  detailed  one  made  by  Suschinskhy  et  al.  [11].  As  in  [1 1],  the 
main  part  of  emission  in  our  study  definitely  belongs  to  the  radon.  In  the  possibility  to 
observe  RB  emission  effect  in  our  case  essential  role  could  be  played  by  higher  sensitivity 
of  our  system  and  observations  of  X-ray  absorption.  One  should  mention  also  a  significant 
difference  in  heights:  the  main  part  of  measurements  in  [11]  were  done  at  the  height  2000 
m,  and  our  —  at  3340  m.  We  emphasize  that  the  mostly  probable  region  of  the  high  values 
of  electric  field  which  can  lead  to  RB  effect  in  thunderstorm  clouds  E~Ec(z)  according  to 
multiple  balloon  measurements  [15,16]  lays  at  the  heights  4-^8  km.  The  characteristic  length 
of  energetic  photon  E  ~  30  h-  100  keV  propagation  in  atmosphere  is  If  <1  km  [6].  So,  the 
probability  to  observe  RB  emission  is  strongly  diminishing  with  height  below  3.5  km. 

We  see  that  though  the  RB  emission  is  effectively  screened  by  radon  it  could  be  brought  to 
light  in  the  ground  based  measurements  as  well.  Of  course  more  detailed  ground  based 
studies  including  higher  heights  and  better  energy  spectrum  and  time  resolution  are  needed. 
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3.6  Conclusion 

In  conclusion  we  will  formulate  briefly  the  main  results  of  our  observations. 

1 .  X-ray  emission  during  thunderstorm  in  mountains  at  the  height  3340  m  is  studied  for  the 
first  time  simultaneously  in  several  different  places. 

2.  It  is  found  that  the  emission  is  well  correlated  at  the  scale  of  the  order  0.5  km  or  more. 

3.  For  the  first  time  well  defined  short  bursts  of  X-ray  emission  are  determined.  The  bursts 
are  highly  correlated  at  the  whole  studied  space  region  (~0.5km). 

4.  The  noticeable  diminishing  of  emission  intensity  penetreating  through  the  13.6  g/cm^ 
carbon  layer  is  established.  It  means  that  the  main  component  of  X-rays  observed  in 
bursts  has  low  energy  (50^80  keV). 

5.  The  X-ray  flux  observed  in  bursts  could  be  estimated  as  10^20  photon/(cm^  sec  sr). 

The  preliminary  analysis  of  our  measurements  allows  to  suppose  that  the  observed 
bursts  of  X-ray  emission  is  the  result  of  bremsstrahlung  amplification  determined  by  the 
runaway  breakdown  effect  in  thundercloud  electric  field. 
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Table  1 


Ground 

Balloon 

Time  duration 

1  5  min 

>  1  min 

Space  correlation 

>  0.5  km 

>  1  km 

Characteristic  energy 

50  -H  80  keV 

60  -  90  keV 

Characteristic  intensity 

10  ^  20  photons/cm^  sec  sr 

40  50  photons/cm^  sec  sr 

Heights 

3340  m 

3500  -  4000  m 
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Figure  3.1 .  Disposition  layout  of  mountain  cosmic  ray  station’s  detectors  (A):  1-  the 
center  of  shower  scincillation  system,  2-  neutron  monitor  and  high-threshould 
gamma  detectors,  3  -  the  system  of  four  WS  low-threshould  gamma-  and  X-ray 
detector  points,  (B):  the  cross-section  of  a  neutron  monitor  unit  together  with  high- 
threshould  and  muon  detectors. 


27 


min-\  1/10*,  min-1  l/10’,  min’  l/10’,  min’ 


^^^^.08.99  19  "U  T. 

LOW  THRESHOLD 
Y-QUANTA 

T— T": 

r 

. 1 . 

- 

1  1  1 

HIGH  THRESHOLD 

- 

“ 

y-QUANTA 

- 

540 

530 

520 

510 

500 

306 

303 

300  W 

297 


li  I  I  I  I  1  1  I  I  I  I  I  I  I  >  I  i . . 1 . i . 1 . 1 . 1 . 1  1 . 1 . 1  i . 1 . t.,  I  I  i  I  t  I  i  I  I  I  ,1  itt  I  I  t  I  I  I  I  I  >  I  I  I 


TIME,  houres 


Figure  3.2.  Typical  behaviour  of  various  cosmic  ray  components  and  atmospheric 
parameters  in  the  absence  of  thunderstorm. 
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Figure  3.3a.  The  behaviour  of  various  cosmic  ray  components  and  atmospheric 
parameters  during  of  the  thunderstorm  at  26  August  1999. 
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Figure  3.3b.  The  behaviour  of  various  cosmic  ray  components  and  atmospheric 
parameters  during  of  the  thunderstorm  at  27  August  1999. 
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Figure  3.3c.  The  behaviour  of  various  cosmic  ray  components  and  atmospheric 
parameters  during  of  the  thunderstorm  at  29  August  1999. 
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Figure  3.3d.  The  behaviour  of  various  cosmic  ray  components  and  atmospheric 
parameters  during  of  the  thunderstorm  at  29  August  1999. 
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Figure  3.4.  Time  dependences  of  correlation  coefficients  Cy  for  three  detector  pairs 
during  the  thunderstorm  on  31  August  1999. 
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where  each  of  i  and  j  indices  designates  one  of  the  low  threshould  detector  points 
(i,  j=1-^4,  nj  and  nj  mean  amounts  of  pulses  obtained  from  bouth  registration 
channels  of  i-th  and  j-th  detector  during  of  1  sec  intensity  counting,  <nj>  and  <nj> 
are  the  corresponding  average  values  defined  over  the  measurements  immediately 
preceding  thunderstorm  moment.  Summation  in  the  above  formulae  was  carried 
out  during  of  100  second  length  time  periods. 
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Figure  3.6.  The  excess  relation  of  low  threshold  to  high  threshold  data 
during  the  thunderstorm  on  31  August  1999. 
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